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I investigate a discrete model of quantum gravity on a causal null-lattice with SL{2, C) structure 
group. The description is geometric and foliates in a causal and physically transparent manner. 
The general observables of this model are constructed from local Lorentz symmetry considerations 
only. For smooth configurations, the local lattice actions reduce to the Hilbert-Palatini action, a 
cosmological term and the three topological terms of dimension four of Pontyagin, Euler and Nieh- 
Yan. Consistency conditions for a topologically hypercubic complex with null 4-simplexes are derived 
and a topological lattice theory that enforces these non-local constraints is constructed. The lattice 
integration measure is derived from an SL{2, C) -invariant integration measure by localization of the 
non-local structure group. This measure is unique up to a density that depends on the local 4-volume. 
It can be expressed in terms of manifestly coordinate invariant geometrical quantities. The density 
provides an invariant regularization of the lattice integration measure that suppresses configurations 
with small local 4-volumes. Amplitudes conditioned on geodesic distances between local observables 
have a physical interpretation and may have a smooth ultraviolet limit. Numerical studies on small 
lattices in the unphysical strong coupling regime of large imaginary cosmological constant suggest 
that this model of triangulated causal manifolds is finite. Two topologically different triangulations 
of space-time are discussed: a single, causally connected universe and a duoverse with two causally 
disjoint connected components. In the duoverse, two hypercubic sublattices are causally disjoint but 
the local curvature depends on fields of both sublattices. This may simulate effects of dark matter 
in the continuum limit. 
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I. INTRODUCTION 

The first order Hilbert-Palatini formulation [1] of classical general relativity (GR) is equivalent to Einstein’s in the 
absence of torsion. Astronomical observations currently cannot distinguish between the two formulations. However, 
in the presence of fermionic matter it seems natural to include an a priori independent SL(2, C) -connection that 
describes the parallel transport of spinors in curved space-time. Including the cosmological term, the Hilbert-Palatini 
action on the four-dimensional Lorentzian manifold M is given by the differential volume form, is [2], 

Shp = ^ f e“ A e" A A - R^^{uj)]eabcd ■ (1) 

Jm 0 

The length scale Ip = here is proportional to the Planck length. Phenomenologically [3] the cosmological 

constant A is positive and incredibly small in natural unit^ with A = AZp = 327rA^|,jg^j^^]^ ^ 1.7 x From a 

semi-classical point of view this small cosmological constant is related to quantum fluctuations of the total 4-volume 
of the currently observable universe (see Appendix]^. The e“ in Eq. 0 are the Einstein-Cartan co-frame 1-forms, 

e“ = eldx\ (2) 

and is the (dimensionless) so(3,1) curvature 2-form, 

+ (3) 

with connection 1-form The Hilbert-Palatini action of Eq.0 does not depend on the frame and is 

defined even if the co-frame is not invertible everywhere. 

This first order formulation in terms of co-frames differs from one in terms of frames in that it is polynomial in 
all fields and depends on the signed invariant volume element. The Lagrangian of Eq. 0 is proportional to det(e) 
rather than |det(e)| and changes sign under improper (local) Lorentz transformations. The action of Eq.Q classically 
is equivalent to the Einstein-Hilbert action only for orientable manifolds with det(e) > 0 everywhere. The Hilbert- 
Palatini action thus includes information on the local orientation of the manifold not provided by the metric. For a 
quantum theory of orientable manifolds, the co-frames of the discretized model must all have the same orientation. 
More important for a discrete version of GR is that forms are geometrical quantities that do not depend on the 
parametrization of the manifold. Contrary to the co-frame e‘^{x) and the connection uj'^{x), the 1-forms e“(x) and 
uj°‘^{x) are geometrical objects. In the first order formulation, the coupling to matter may also be expressed in terms 
of forms only (see Appendix . The remaining local internal symmetry of such an oriented model is 5'0(3,1) and 
causality further restricts this to the connected component 50+(3,1) ~5L(2, C). 

In the absence of torsion, the Hilbert-Palatini action of Eq.Q is (twice) the real part of the action for self-dual 
connections and the classical phase space may be restricted accordingly [3HS]- This has been the starting point of most 
Hamiltonian approaches to quantum gravity[ll [SHZ]. 

It is difficult to regulate Hamiltonian continuum theories. A lattice formulation in terms of geometrical quanti¬ 
ties potentially could overcome these difficulties. However, in dynamic triangulations based on Regge’s simplicial 
description [SI in] of GR, a crumpling instability leads to non-causal conhgurations [TDlll3j . Ambjprn, Jurkiewicz and 
Loiinn nni realized that a restriction to causal manifolds may stabilize this model in some regions of the coupling 
space. They solved the causality constraints |lbl ITT] with a particular dynamic triangulation. 

Here I elaborate on a recently proposed [T8| causal lattice formulation of Eq.Q. The discrete counterparts of 
forms provide a geometric triangulation of space-time. To accommodate the conventional coupling to matter, the 
co-frame 1-form will be treated as distinct from the SL{2, C) connection. The co-frame thus will not be considered 
a component of a single gauge connection as in [T3H33] . 

The resulting model is a hybrid between a conventional lattice gauge theory and the rigid simplicial approach of 
ReggejSl in], Ambjprn and Loll[T31 |TH HSl |23]. It differs from ordinary lattice gauge theory in that the geometry of 
the lattice depends on the field conhguration on it, but unlike the Ambj0rn-Loll approach, the basic variables are 
5L(2, C)-variant spinors. The lattice retains a fixed (hypercubic) coordination but its cells vary in shape and size. 
Distances between nodes of the lattice in particular depend on the configuration of lattice variables. We consider 
triangulations for which the separations between neighboring events are light-like. Such a null-lattice does not of itself 
introduce a length scale and an ultraviolet cutoff is required to ensure a minimal coarseness of the triangulation. The 


^ In natural units h = c = Ip = 1. The Minkowski metric r) = diag(l,..., 1, —1) with ’’mostly positive” signature (-1- -|- H—) is used 
throughout and Einstein’s summation convention for repeated diagonal indices is adopted. The symbol £{abcd) = £abcd = gives 

the sign of the permutation of its arguments with e(1234) = 1. Following common conventions, lower- (upper-) case letters from the 
beginning of the Latin alphabet denote Lie-algebra indices of so(3,1) and sl{2,C) respectively. Latin indices from the middle of the 
alphabet are used for tensors. Greek indices will be reserved for labeling parametrization invariant (lattice) objects. 
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latter breaks scale invariance dynamically but allows one to compute finite distances on a finite lattice. Measured 
quantities are related to lattice amplitudes with certain geometrical characteristics. The critical limit of this model 
thus differs from that of ordinary lattice gauge theories where the lattice geometry is assumed fixed from the outset. 

The present model is based on the same principle as the Global Positioning System [251 - 127] (GPS): The intersection 
of forward light-cones from four spatially separate events is a later event. 

In the GPS, the (spatial) space-time separations of the emitted signals are determined using atomic clocks on 
satellites with known orbits. These satellites furthermore emit a string of signals (events). In the idealized lattice 
setting, the emitting event is itself specified by the intersection of four previously emitted signals. The intersection 
of light cones in this sense defines an event of the causal null lattice. The resulting network of events with light-like 
separation in principle could be traced all the way back to the Big Bang, the single event from which all others derive 
causally. 

The analogy with the GPS implies that exactly 4 events must lie on the backward light cone of every event and that 
every event must illuminate 4 others. The coordination of the network thus has a fixed value of 2d in d space-time 
dimensions and is not arbitrary. We will use the nodes of a hypercubic lattice to label the intersection events and 
indicate their causal relation. 

The use of light-like links solves the causality problem. However, instead of reality constraints|Sl|B], certain condi¬ 
tions have to be met to ensure that a configuration gives a triangulated manifold. These consistency conditions are 
derived in Sect|^ and a local Topological Lattice Theory (TLT) which enforces them is constructed in Sect, 

The model thus represents causal manifolds by simplicial complexes. Its point set of events preserves causal order, a 
feature this model shares with causal set theorv ^BIIdOj . However, local Lorentz invariance is imposed by averaging 
over equivalent configurations rather than emergent and the network is designed to represent triangulated causal 
manifolds only. In the strong coupling limit of Sect |IX| the model becomes a random measure of spatial lengths that 
satisfy certain constraints. 


VB 


II. THE LATTICE 


We will consider a four-dimensional topological lattice A with hypercubic coordination. This lattice labels events 
and gives causal relations between them, but does not of itself specify a manifold. Topologically hypercubic implies 
that the adjacencies of events are those of a hypercubic lattice but, of itself, does not specify distances or angles 
between neighboring events. It does imply that each event has 8 neighbors and that links between neighboring events 
can be oriented. 



FIG. 1. a) [Wikipedia File:HC Pl-P3.png] Tetrahedral-octahedral triangulation of a spatial hyperplane. Each node has 
coordination 12. b) the immediate neighborhood of a node on the hyperplane. The four B-tetrahedrons with a common node 
are shown in blue and the four F-tetrahedrons are shown in yellow. Each of these two sets of tetrahedrons includes all nodes 
of the spatial hyperplane although only their union with the interstitial octahedrons is space-filling. Every edge is shared by 
one F- and one B-tetrahedron and every vertex is common to four B- and four F- tetrahedrons. Note that the 8 faces of the 
octahedral voids coincide with those of four B- and four F-tetrahedrons. For an general 3-dimensional Cauchy surface the 
tetrahedrons in general will be neither equal nor regular. 
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By analogy with the GPS, the coordination of the null-lattice in 3-1-1 dimensions is 8 , with 4 future and 4 past 
events associated with each node. This may be realized by a lattice with hypercubic structure. That this is the only 
possibility can be seen by examining flat Minkowski space-time. To this end consider first a spatial hyperplane Pq of 
Minkowski space. The intersection of light-cones from a tetrahedron of adjacent nodes in Pq defines a single node in 
a later hyperplane. To conserve the number of nodes, each node of a spatial hyperplane therefore must be common 
to 4, otherwise disjoint tetrahedrons. One thus is led to the triangulation of Pq by a regular tetrahedral-octahedral 
honeycomb. The nodes of the hyperplane Pq (at ’’time” 714 = 0([^thus can be labeled by even integers, 

Po = {(2ni, 2n2,2n3,0); ni, 712,713 G Z with ni-I-712-I-713 = 0 mod 2} . (4) 

This tetrahedral-octahedral tesselation of flat 3-dimensional space with coordination 12 is shown in Fig. The 12 
spatial neighbors of a node form 4 disjoint sets of 3 that neighbor each other. Together with the common node, these 
are the vertices of 4 (spatial) tetrahedrons. We will call this a set of F-tetrahedrons (or ’’forward” tetrahedrons). 
In Fig. [^) they are colored yellow. The complementary set of 4 (blue in Fig. [^)) B-tetrahedrons (or ’’backward” 
tetrahedrons) have the same common node and every edge of a B-tetrahedron also is the edge of an F-tetrahedron, 
but B- and F- tetrahedrons do not share a face. The vicinity of the node (0, 0,0,0) of the lattice A includes the 
following sets of nodes. 


(-1, 1, 1,-1) (0, 0, 0. 

( 1 .- 1 , 1 ,- 1 ) illumina.tion ^ (2,-2, 0, 

( 1, 1, —1,—1) intersection (2, 0,-2, 

(-1,-1,-1,-1) (0,-2,-2, 


(i.-i.-i, 


B-tetrahedrons 


0 ) (- 2 , 2 , 0 , 0 ) (- 2 , 0 , 2 , 0 ) ( 0 , 2 , 

0 ) ( 0 , 0 , 0 , 0 ) ( 0 ,- 2 , 2 , 0 ) ( 2 , 0 , 

0 ) ( 0 , 2 ,- 2 , 0 ) ( 0 , 0 , 0 , 0 ) ( 2 , 2 , 

0 ) (- 2 , 0 ,- 2 , 0 ) (- 2 ,- 2 , 0 , 0 ) ( 0 , 0 , 

illumination fj, intersection 

1 ) (- 1 , 1 ,- 1 , 1 ) (- 1 ,- 1 , 1 , 1 ) ( 1 . 1 , 


2. 0) G 

p’ pj > F-tetrahedrons 

0, 0) J 

1. 1) . (5) 


If the four nodes of a column/row label the vertices of a B/F- tetrahedron, Eq.([^ may be interpreted as follows: 
the forward light cone^of the B-tetrahedron at depth 714 = —1 illuminate the four F-tetrahedrons (rows) at depth 
774 = 0. The events of the selected B-tetrahedron of the 714 = —1 hyperplane thus are on the backward light cone of 
the (0, 0,0, 0) node. Each node at 774 = —1 is at the intersection of four backward light cones based at the vertices of 
an F-tetrahedron in the 774 = 0 hyperplane. The rows/columns of the matrix of vertices on the 774 = 0 hyperplane 
label the vertices of F/B- tetrahedrons with the common node (0, 0,0,0). The forward light cones from the vertices 
of a B-tetrahedron in the 774 = 0 hyperplane intersect at a vertex of the 774 = 1 hyperplane and the corresponding 
B-tetrahedron is illuminated by this vertex. The four vertices obtained as intersections of light cones from four B- 
tetrahedrons that include the (0,0, 0, 0) node form the F-tetrahedron of the 774 = 1 hyperplane in Eq.([^. Its vertices 
are on the forward light cone of ( 0 , 0 , 0 , 0 ). 

Although but a small local section of the whole lattice, Eq.([^ illustrates the general construction: 


• A spatial Cauchy surface is triangulated by B-tetrahedrons with no common edge. Every vertex of this trian¬ 
gulation is common to four B-tetrahedrons. 


• Vertices of the next spatial hypersurface correspond to the intersection of four forward light cones emanating 
from the nodes of a B-tetrahedron of the earlier Cauchy surface. 


• The vertices obtained from four B-tetrahedrons with a common node form an F-tetrahedron of the later spatial 
hypersurface. They lie on the forward light cone of the node that is common to the B-tetrahedrons. Four of 
these F-tetrahedrons have one node in common but do not share an edge. 


• The corresponding B-tetrahedrons of the later hyperplane share an edge with each of six surrounding F- 
tetrahedrons. 


• the procedure is repeated to obtain the spatial hypersurface of the next depth. One thus arrives at a foliated 
causal description of space-time. 

Fig. [2] illustrates the analogous construction in 1+1 and 2+1 space-time dimensions. The hypercubes in 2+1 
dimensions are three-dimensional cubes and the tetrahedrons of a spatial hypersurface become triangles of a two- 
dimensional, hexagonally triangulated surface. Note that the F- and B- triangles individually suffice to determine the 


^ In Minkowski space the spatial hyperplanes correspond to global time, but a spatial hypersurface with constant depth 714 of Lorentzian 
space-time in general will not be one to constant cosmic time. 

^ To represent Minkowski space by an integer lattice the speed of light may conveniently be taken as c = l/y/S. 
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(hexagonally) triangulated two-dimensional surface. In 1-1-1 dimensional space-time the cubes reduce to squares and 
the (F/B) triangles become line segments that coincide. Note that the number of independent spatial lengths is the 
number of geometrical degrees of freedom in any dimension, that is 1,3 and 6 per node in d = 2, 3 and 4 dimensions 
respectively. 




a) b) 

FIG. 2. Topological lattices in 1 -I- 1 and 2-1-1 space-time dimensions. Forward null vectors are in red and the corresponding 
backward (B) and forward (F) spatial simplexes are shown in blue and yellow respectively, a) In 1-1-1 dimensions these spatial 
simplexes are 1-dimensional and coincide. Forward and backward light-cones are triangles with indicated apexes. Events and 
their causal relation correspond to nodes and links on this topologically square null-lattice, b) In 2-|-1 space-time dimensions the 
topological null-lattice is cubic and the 2-dimensional spatial hypersurfaces are hexagonally triangulated with coordination 6. 
[As discussed in the text, the construction of a null-lattice in 3 -I- 1 dimensions is analogous, with hypercubes replacing cubes 
and tetrahedrons replacing the triangles. The tetrahedral-octahedral triangulation of 3-dimensional spatial hypersurfaces in 
Fig. [^has coordination 12.] The 3(4) vertices of blue B-triangles (tetrahedrons) are the apexes of future light cones (with rays 
shown in red). Their intersection corresponds to the node of a yellow F-triangle (tetrahedron) on a future Cauchy surface. The 
subset of F- or B-triangles (tetrahedrons) both triangulate a 2(3)-dimensional Cauchy surface. The triangles (tetrahedrons) of 
a causal 2(3)-1-1 dimensional Lorentzian manifold in general are not regular and the lattice is deformed as in a). 

We triangulated flat Minkowski space, but the construction can be generalized to causal Lorentzian manifolds. 
Future events again are at the intersection of (four) forward light-cones from a spatial B-tetrahedron of earlier events, 
although these B-tetrahedrons generally are not regular. Four spatial B-tetrahedrons with a common event define 
the events of a future F-tetrahedron. The construction fails if the four light cones emanating from the vertices of a 
B-tetrahedron do not intersect. With only a finite number of nodes this can occur near curvature singularities such 
as the singular world line of a black hole or the singular surface of a cosmic string. In the present construction such 
singular defects can only occur in the critical limit where the number of nodes becomes arbitrary large while the total 
4-volume remains finite. 

The main challenge will be to ensure that field configurations describe causally triangulated manifolds. Since each 
space-like edge of a B-tetrahedron is shared with an F-tetrahedron, the endpoints of every spatial edge lie on the 
backward light cone of a unique future event and on the forward light cone of a unique past event. One can glue the 
edges of F- and B-tetrahedons to a pure 3-complex only if these (spatial) lengths coincide for the respective past and 
future 4-simplexes. The Topological Lattice Theories (TLT’s) constructed in Sect|v| and Appendix [P] ensure this 
consistency condition is satisfied and allow one to uniquely construct the triangulated causal manifold. 

Eq.([^ implies that events of the piecewise linear oriented manifold may be labeled by integer quartets of a 
hypercubic lattice, 


4 

A = {n = ^ , 


( 6 ) 


with forward lattice displacements. 


Ai = (I,-l,-I,l), A2 = (-!,!,-1,1), As = (-l,-1,1,1), A4 = (1,1,1,1) . 


(7) 










































6 


The labels of 12 spatially neighboring events of a node on the same hyperplane are found by adding the (12) differences 
of these displacements to the label of the node. This completes the construction of the topologically hypercubic lattice 
A that labels events of the causal universe. 

In an alternate scenario connections and co-frames are placed on separate links. In addition to A we in this case 
consider the lattice with ’’time-conjugate” forward displacements, 

Ai = (1,-1,-1,-1), A2 = (-i,i,-i,-i), A3 = (-i,-i,i,-i), A4 = (i,i,i,-i) . (8) 


By construction these satisfy. 


- Ai, -I- Ai, - = 0 for all v 

A^ -I- Aj, -I- Ap -I- Act = 0 for {/X, i/, p, tj} = { 1 , 2 ,3,4} . 


(9) 


The displacements of Eq.(|^ and Eq.Q are oppositely orientecj^ and related by the orthogonal matrix O with 
det O = —1, 


( 10 ) 


1 

1 — 1 

1 

-1 

-1 

1 — 1 

1 

1 

-1 

-1 

1 — 1 

1 

1 — 1 

1 

1 

-1 

1 — 1 

1 

1 — 1 

1 

-1 

1 


The conjugate forward displacements of Eq.(|^ generate a hypercubic lattice with nodes labeled by, 

4 


A = {n = ^ npAp; S 


( 11 ) 




By Eq.§, the labels of all ’even’ nodes of A and A coincide and the corresponding events are identified, with 

4 4 

A n A = {n = ^ UpAp; ^ Up = 0 mod 2} . (12) 






However, the union AU A leaves ’odd’ nodes of A out in the cold and generally leads to an unstable model. This can 
be avoided by including a second, translated copy h! of A that coincides with A at odd nodes, 

4 

A' = {n = A 4 -H ^ ripAp; rip e Z} . (13) 

/i=i 


The combined lattice. 


Ad = A U A' U A = A U A' , (14) 

is the union of the hypercubic lattice A with vertices of the A' lattice at the centers of its cells. This body-centered 
hypercubic Cl lattice, the 16-cell honeycomb of 4-dimensional Euclidean space, is the 4-dimensional analogue of a 
body-centered cubic lattice. The Euclidean distance to each of the 24 nearest neighbors of a vertex on this lattice 
is 2 in the present parameterization and unit 3-spheres whose centers coincide with vertices of A^ have the densest 
packing in 4-dimensional Euclidean space[3T]. The hypercubic sub-lattices A and A' are in dual positions, i.e. A 
could be interpreted as the dual hypercubic lattice to A'. A C Ad is a hypercubic sub-lattice of Ad whose even and 
odd nodes coincide with nodes of A and A' respectively. 

Inclusion of the dual hypercubic lattice leads to a more symmetric and appealing construction at the price of 
’’doubling” the number of causally connected components. Locally, the two components, A and A', of the ’’duoverse” 
Ad are causally disjoint although the space-time geometry of each connected component will be determined by the 
field content on both. From the point of view of one of its causally connected components, the duoverse Ad resembles 
a single universe constructed on A only. The (gravitational) field content of both scenarios in fact is the same. 


^ The lattice displacements {Ap} ({Ap}) form left(right)-handed systems. 
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III. FIELDS AND OBSERVABLES 


To obtain the Lorentzian lattice model with causal dynamics, consider first the topologically hypercubic lattice A 
of the universe. For the duoverse a similar construction holds on the hypercubic A sub-lattice. 

Links are ordered pairs of neighboring nodes of A. We adopt standard conventions and use [n, n'] or [n, n + A^] = 
[n, n + fi] = [n,/r] to identify links of A - whichever notation is more convenient. A is naturally oriented by the 
displacements of Eq.Q. Links [n -I- /i,n] = [n -|- are said to be reversed. Sometimes more general paths 

[no, rii, n 2 ,... ] have to be considered. They are given by a contiguous set of oriented links [n^, n^+i] e A. 

The variables (n) of this lattice formulation are a finite lattice version of the continuum I-form {n)dx^ 
in spinorial representation [5^ . The differential dx^ is represented by a contravariant vector that approximates the 
’’displacement” to the node n + g, along the geodesic C[n, n -|- g\ between the two events. 


E^^\n) := {n)E{n) ^ e 


AA', 


AA 


,dx^ 


Ip 


(15) 


The contravariant displacement vector is uniquely defined using the affine parametrization of a geodesic with 
a;(A = 0) = a:(n) and a::(A = I) = x{n + g) and AA = 1. In natural units it is the tangent u^(A = 0) = u{((n) to the 
so affine parameterized geodesic C[n, n -|- fi] joining the two event^ 



1 dx'^ 
Ip dX 


(A = 0) 

C[n,n+/i] 



(16) 


The dimension of the displacement vector has here been absorbed by so that all variables of the lattice model are 
dimensionless. We effectively are measuring length in units oi Ip = \J’yi'xhG j « 1.6 x time in units oilpjc 

and energy in units of hc/lp. Ip here is a dimensional unit on the same footing as h and c. The only gravitational 
coupling of this lattice model is the dimensionless cosmological constant A = AZp. It plays a role analogous to that 
of the (in 3-1-1 dimensional space-time) dimensionless couplings a = j(he) and g^ of the electro-weak and strong 
interactions. 


The displacement vector lf)(n) of Eq.(16) for given events n and n-|-^ is a (uniquely defined) contravariant vector 


and its contraction with the co-frame ef^(n) to E„(n) in Eq.(15l is a scalar. A smooth change of coordinates in 




the vicinity of the event n in particular transforms l^(n) and the co-frame efc(n), but not E^('a). For lack of a better 
name, we refer to {i?p(n),/r = 1,... ,4} as a lat-frame at the node n, even though this ’’frame” is parametrization 
invariant. 

In Appendix first order actions for spinors, scalars and gauge fields are written in terms of (diffeomorphism 
invariant) forms. Displacements and co-frames thus never appear separately and the corresponding actions for the 
discretized model can be written in terms of lat-frames only. This is desirable, since the diffeomorphism group 
does not act on the triangulation of a continuous manifold by a finite number of points. Apart from an internal 
SL(2, C) symmetry, the present lattice model thus is a dynamical triangulation that preserves coordination number. 

The matrices Efj^(n) transform homogeneously under local S'L(2, C) transformations. 


Ef,{n) —^ 5(n)F;^(n)g1'(n) with ^(n) e SL{2,C) , 
and the SL(2, C) invariant inner product on this space of anti-hermitian matrices is. 


(17) 


X-Y ■= -iTieX'^eY = - sabY^^'ea'B' 


e R 


(18) 


where = —X* denotes the transpose of the anti-hermitian matrix X = —X^ and e = (J) is the SL{2, C) invariant 
tensor. 

It sometimes is more convenient to consider lat-frames as so“'"(3,1) vectors. The relation to their spinorial repre¬ 
sentation is provided by an orthogonal basis of anti-hermitian 2x2 matrice^ 

{erf^ ,a = I,...,4; tJa = -cr^ with Tva^ab = 26^, where a a>b := SBCC^a^ sc A'}- (19) 


^ This geometric interpretation of (n) is owed[33] to T. Jacobson. 

® We in particular consider a basis in which the a = (cri, C72, cr 3 )-matrices are traceless anti-hermitian generators of an su{2) algebra and 
C74 = 2l. 








The spinor and vector components of lat-frames are related by, 




(20) 


and the scalar product of Eq.(18l in the vector representation is, 

X • r = X^TIabX'^ 


( 21 ) 


with ? 7 ab = diag(l, 1,1,-1) . 

As discussed previously, to ensure causality and in analogy with the GPS, all forward displacements are chosen 
light-like. 


-£^^i(n) • E^(n) = gjfc(n)P (n)l^(n) = 0 with -iTvE^>Q. 


( 22 ) 


The same causal manifold may be triangulated in a number of ways, but its triangulation by a topological null- 
lattice is unique up to the placement of events on the boundary of the manifold. [For a more detailed discussion of a 
conic manifold see Appendix]^] 

Up to S'L(2, C) transformations, a null lat-frame represents the 6 spatial distances > 0 between four events 

on the forward light-cone of n. 


(n) = - E,) ■ {E^ - E,) = -2E^{n) ■ E,{n) = -2g,fc(n)Z;(n);^(n) > 0 . 


(23) 


These lengths are coordinate invariant geometrical quantities. Modulo SL(2, C) transformation, a null lat-frame thus 
describes 6 geometric degrees of freedom. 

In general, lat-frames are not proportional to the co-frames in any coordinate systenj^ However, they may be 
identified with (forward) null co-frames in Minkowski space-tim^ for instance, 


i A-.)' = ■ ( 


V3-1 1+i A 

1-i V3-I-1 J ’ 


= I 


■ f V3-1 -1-i 

V -1+i V3+1 


(24) 


5'L(2, C)-transformed and scaled null lat-frames give an equivalent description of discretized Minkowski space-time. 
Note that the lattice displacement vectors in S' are l^ = aS^. The global lattice constant a gives the coarseness 
of the discretization (in units of Ip). This reflects the possibility of choosing a local inertial system that describes 
the immediate neighborhood of an event. The resulting piecewise linear approximation to a smooth manifold can be 
justified for a sufficiently ’’fine” lattice, but fails to reproduce curvature singularities. These singularities may only be 
recovered in the critical limit where the number of events becomes infinite. 

Eq.(18) and Eq.(22) imply that null lat-frames are singular matrices and may be represented by complex bosonic 
2-component spinors {A = 1,2),- 


® C)n" = (n) 


f * \AB' 


f-*B' I 


(25) 


where denotes the complex conjugate spinor. A primed upper case Latin index indicates that the spinor transforms 
with the complex conjugate representation of SL{2, C). There is no summation over the repeated Greek inde^l^/i in 
Eq.(25|. Gomponents of the conjugate y'’ of any spinor y are. 


Xa '■= X £ba ■ 

^^(n) and its conjugate ^()(n) thus transform inversely under SL{2, C), 

^' 7 l 3 (n)C;f(n) ^ (n), for ^(n) S S'L(2, C) . 

A lat-frame in addition is invariant under local [/“^(l) phase transformations of the spinors, 

^^(n) ^ ^(,(n) = and ^*{n) -b ^'*{n) = (n), V’/x(n) S R . 


(26) 


(27) 


(28) 


^ The displacement matrices ^^(n) generally do not correspond to a smooth coordinate transformation dx^/dx'^{n) at more than one 
node. 

® The null lat-frames for instance can be taken proportional to co-frames of the system S' related to Minkowski space with metric 
= diag(l, 1,1, —1) by the linear coordinate transformation 4x'^ = —x^ x"^ x^ = x^ — x"^ x^ x'^/y/^ Ax'^ = 

x^ + x‘^ — x^ + fc‘^/\/3, Ax"^ = —x^ — x^ — x^ + x‘^l\fZ. This transformation preserves orientation and, as required by Eq.( |23| l, scalar 
products of the lat-frames are negative semi-definite. The metric in S' is ^ ^ y and vanishes for fi = ia. 

^ Only diagonally related, repeated indices, i.e. ^re summed. 
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A spinor may be compared to another by parallel transport along links of the lattice to a common node. The 
parallel transport of a spinor from n + to n is given by matrix C/^(n) = t7[n, = C/[n, n + A^] € SL{2, C). Nodes 

of the lattice thus are associated with four spinors and its oriented links [n, /i] with parallel transport matrices. On 
a (periodic) hypercubic null-lattice with N nodes there are altogether AN spinors and AN transport matrices. Under 
the SL{2, C) structure group the transport matrices transform as (^(n) G SL{2, C)), 


^ 5l;(n)U^r,(n)5 ^ ^(n -t- fi) . 


(29) 


Although consistent with Eq.(29), the reversed lattice link will not be associated with the inverse-, but rather with 
the hermitian conjugate transport matrix, 


C/[n, n'] = n] . 


(30) 


Since the representation of SL{2,C) is not unitary, Eq.([30| eliminates closed (Wilson) loops of S'L(2, C) transport 
matricefrom the S'L(2, C)-invariant observables one can construct for this lattice. Eqs. (27), (29) and (30) imply 
that the basic SL{2, C) invariants of this model are contiguous strings of transport matrices bookended by spinors, 

4^^[no,ni,...,nr] := ^))(no)C/[no, ni]C/[ni, ns]... 17[n^_i, n^]^i.(n,.) , (31) 

where [ng, ni,..., n^j is a continuous chain of forward oriented links. For reference we separately list and name the 
three shortest invariants. 


X/ii/po-(n) 


:= =: 4°^[n] , 

:= CA(n)C//B(n)Cf (n-f p) = C;i(n)C/[n, ni]C,,(ni) =: /W[n,n-kp] , 
:= ^M(n)C//B(n)C/J^c(n + (n + p + a) 

= Cp(n)U[n,ni][7[ni,n2]C^(n2) =: ll^J[n,n + p,n + p + a] . 


(32a) 

(32b) 

(32c) 


The invariants of Eqs. (31) and (32) in general are complex and physical observables necessarily depend on the 
corresponding complex conjugate invariants as well. The path of the latter is reversed and the y a re bookended 
by hermitian conjugate spinors. Invariance under the local U^{1) phase transformations of Eq.(28) implies that 
phy sica l obser vable s locally conserve four separate spinor numbers. They thus depend only on the null lat-frames of 
Eq.(25|. Sect V B| exploits the symmetry to construct a local TLT that constrains the space of lattice spinor 

configurations to those that correspond to triangulated causal manifolds. 

The magnitude of the symplectic form /pi/(n) = —/,,^(n) on the space of spinors at site n is the spatial length 
el,{n) of Eq.@, 


^L(n) = TreEUn)eE^{n) = |/p,.(n)p > 0 , 


(33) 


where we used Eqs. (25) and (32a). Note that the length vanishes only if the two spinors and are linearly 
dependent. 

In four dimensions the 6 complex variables of the anti-symmetric matrix are constrained by the fact that 

its Pfaffian vanishes. 


Pf(/(n)) ^ XI £(M^P0-)/M>^(n)/pc7(n) = /i2(n)/34(n) -k /i3(n)/42(n) -k /i4(n)/23(n) = 0 . 


(34) 


The matrix /(n) thus depends on 10 real parameters only. Four of these are overall phases of the spinors [32] (see 
Eq.(|^). The remaining 6 real degrees of freedom are the spatial lengths fpi/(n) > 0 of Eq.(33). 

The vanishing Pfaffian implies that the magnitudes of the three complex numbers in Eq.(34) satisfy triangle 
inequalities. 


0 < a(n) < 5(n) -k c(n), 0 < 6(n) < c(n) -k a(n), 0 < c(n) < a(n) -k 5(n) , (35) 


The identification f7[n,n'] = U ^[n',n] would include closed (Wilson) loops of the form Trl7[ni, n 2 ]P[n 2 , na]... 17[nr, ni] as basic 
SL(2, C) invariants. 
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with, 


a(n) 

6(n) 

c(n) 


=£i2(n)£34(n) = |/i2(n)/34(n) 
=£i3(n)£24(n) = |/i3(n)/42(n) 
=£i4(n)£23(n) = |/i4(n)/23(n) 


For non-negative a, b, and c, the triangle inequalities of Eq.(35) can be combined to the single inequality, 


0> a'^+ + c‘^ - 2a^b^ - 2b^c^ - 20 ^ 0^1 = detil^in) . 

— I n ' 


(36) 


(37) 


The triangle inequalities of Eq.(35) thus are equivalent to requiring that the local four volume E(n) = det ^^^(n) 

be real. The fact that causality restricts the spatial lengths £^jy(n) so simply is one reason for considering null lat- 
frames. The converse holds as well: given six non-vanishing positive lengths >0, l</i<z^<4} with 

det^^y < 0, an anti-symmetric complex matrix with = ijiv can be constructed whose Pfafhan vanishes. The 
spinors then are reconstructed as follows. 


The vanishing Pfafhan of Eq.(34) implies that /^iy(n) and its dual /^“^(n) := antisym¬ 


metric 4x4 matrices of rank 2. The null space of /^*^(n) thus is spanned by two linearly independent vectors 

S/j, ■) 




,nn)C^in) = 0, for A = 1,2 


(38) 


Up to normalization these dehne (n) itself. The fact that any 4-dimensional anti-symmetric matrix with vanishing 
Pfafhan is represented by a set of spinors is exploited in Sect|V|and Appendix [Cl 

IV. FIELD ASSIGNMENTS AND LATTICE ACTIONS 

The model is completed by assigning this held content to links and nodes of the lattice and constructing the lattice 
action and integration measure. We will consider two distinct scenarios: 

A. The Universe 


Only the hypercubic A lattice of Eq.Q is utilized in this case. Each of its oriented links [n, fj] is associated with 
an SL(2, C) tran spor t matrix Up(n). The reversed link is associated with the hermitian conjugate transport matrix 
[/^(n) as in Eq.(30). 

Each node n furthermore is associated with a null lat-frame {i?^(n), /i = 1,... ,4} that gives the light-like 
displacements in the forward direction from n to the events {n -f A^, ^ = 1,..., 4}. 

The lattice action is constructed from purely imaginary SL(2, C) invariant lattice 4-forms with real coefficients. The 
two most local candidates for the lattice action are, 


Vpiypain) ■—TiEfj,{n)eE^{n)eEp{n)eEj{n)e = 2ie{fiiypa) det[E] := 2is{pi^pa)V{n) 
=fpA^)Upi^)fpai^)f<yp{^) , 

Vp^^pain) :=xpupain)x*p^^pin) = TreUj(n)eUp(n)C/<,(n -k p)E^{n + p + cr)C/^(n -k cr)C/](n) 

=CA(n)U/B(n)Ufc(n + p)^^{n + p + a)ff {n + a + p)C/*/(n -k a)U*^{n)(,*^P{n) 


These contributions to the action are depicted graphically in Eigs.[^) and[^) respectively. 

The other two quasi-local observables, 

:=Tre£;J(n)eE^(n) = |/^^(n)|^ = elp{n) 

■ IXfiiypa- (n)p = TTeE^{n)sUp{n)U^{n + p)E^{n + p + cr)t7](n -k p)U^(n) 
=CM(n)C^/B(n)Ufc(n -k pA^ (n -k p -k (j)C^{n -k cr -k p)U*^{n + p)U*/{n)C^\{n) 


(39a) 

(39b) 

(40a) 

(40b) 


are real. The first of these invariants are the spatial lengths of Eq.(33) and is symmetric in its indices. The totally 
anti-symmetric part of Q^ypcr(n) is the lattice analog of a topological contribution to the continuum action. This 
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FIG. 3. Quasi-local contributions to the invariant action on a null-lattice in d = 3 -I- 1 dimensions, a) The cosmological term 
Vui/po-(n) at site n of Eq.(39al. Spinors are shown as red segments alo ng l attice direction /r, hermitian conjugate spinors 

are depicted in green. Color flows from conjugate spinors dehned by Eq.( 261 to spinors and from to . Two directed 
line segments of the same color joined at n form the SL(2, C)-invariant (red) /pi/(n) of Eq.(32al (or its complex conjugate 
(green)). The contribution to the cosmological term of site n is the product of the se invariants with the sign determined by 
the permutation of {^upa). b) The Ppiypa-(n) curvature contribution to Eq.(39bl. The depiction of spinors at sites n and 
n -f p -I- (T is as in a), but these now are parallel transported along links in p- and a- direction by SL{2, C) transport matrices as 
shown. The (red) spinors are parallel transported by the SL{2, C) -matrices U, whereas their (green) hermitian conjugates 
counterparts are transported by the hermitian conjugate S'L(2, C)-matrices U~^. The contribution to the lattice 

Hilbert-Palatini action is given by the product of these two overlaps with a sign determined by the permutation of {pvpa). c) 
The Q^vpa topological term of Eq.(40bl. In the continuum limit his term does not affect classical equations of motion and 
could be included in the lattice action with a purely imaginary coefficient. Similar to the cosmological contribution a), the 
corresponding path does not enclose a plaquette and does not measure local curvature. 


term could be (and its continuum analog often is[4l |34l [35]) included but will not be considered in the following. 
We furthermore do not investigate other, less local, contributions to the lattice action. The invariant (real) lattice 
actiorp] for a causally connected universe thus is, 

•^HP' = * 51 51 + ^''^M^P'r(n)] , (41) 

^yp<7 n^A 

where A is the dimensionless cosmological constant. We have absorbed the coefficient (and dimensionality) of the 
curvature term in the normalization of the spinors (or equivalently, of the lat-frame E^). On a hypercubic lattice 
the HP-action takes the same form whether the lat-frame is light-like or not, but the relation between events of A in 
general would be acausal. Null lat-frames assure causality, but only a subset of such lattice configurations corresponds 
to triangulated manifolds (see SectQ. 


B. The Duoverse 


The duoverse is obtained by placing the held content on the lattice of Eq. (14). In some ways this scenario 
is more restrictive and appealing than the universe. In analogy with optics, nodes n € A of the hypercubic sub¬ 
lattice A C A.d will be referred to as ’’active” nodes. All other nodes h G Ad\A of A^ are ’’passive” or inert in this 
construction. In Table ||| and Fig. active nodes are denoted by solid dots and passive ones by open circles. Half 
the nodes of A and A' are passive. Active links are between two active nodes whereas passive links are between an 
active and a passive node. The hypercubic sub-lattices A and A' are composed of passive links only. Depending on 
whether the initial or hnal node is active, oriented passive links are of two kinds [n, h] and [h, n] . There are no links 
between two passive nodes in this construction. 5'L(2, C) transport matrices Up^ijv) reside on oriented active links 
[n, n' = n -P A ^] S A only. The reversed active link being associated with the hermitian conjugate transport matrix 
as in Eq.(30) and TableFig. [^depicts these assignments for a pcr-plaquette of A (or A') and a /^/-plaquette of A 
that have two active nodes in common. 

The number of degrees of freedom of the universe and of the duoverse are the same, since they have the same number 
of active nodes and active links. In both arrangements each active link is associated with an SL{2, C) transport matrix 


The weight e®’®HP of a configuration is a pure phase. 
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oriented 

reversed 

active links 

Ap , 

n •-• n 

t/e(n) 

C^^(n) 

passive links 

n •- y o n 

Cp(n) 

^i(n) 

Ap 

n o-^ • n 

ep(n) 

C?(n) 


TABLE I. Links and field assignments of the duoverse. Solid dots (•) denote active nodes of the sublattice A, whereas open 
circles (o) are passive nodes of A (or A'). Assignments depend on whether a closed path traverses a link in the direction of 
its orientation or reverse to it. Passive (o) nodes act like mirrors: they are never crossed and reverse the direction of a path. 
A path on a passive link thus always retraces itself. Active nodes on the other hand do not reverse the direction of a path. 
Physical observables of the duoverse are associated with closed lattice loops. 



FIG. 4. Plaquettes of the duoverse. A ^i/-plaquette of events (solid nodes) in A and the corresponding pa-plaquette of A (A'). 
Solid dots (•) denote active events of A. Passive events of A (A') are shown as empty circles (o). Note that spinors ^ (^‘^) 
reside on outgoing (incoming) passive links whereas SL(2, C) transport matrices occupy active links. Links traversed in reverse 
direction are associated with hermitian conjugate variables. The local SL(2, C) structure group acts at active sites of A only. 


and each active node with four forward null vectors. The structure group SL{2, C) acts at active nodes only. The 
passive nodes of the duoverse carry no additional dynamical information. In the next section the triangulated manifold 
is reconstructed from outgoing passive links of A (A') only. In the duoverse forward and backward light cones of 
active nodes are glued at passive nodes, whereas they are glued at active nodes in the universe. The universe and 
duoverse thus differ only in the manner the complex is formed from the basic simplexes. 

A path on A either traverses a link in the direction of its orientation or in the reverse direction. In the duovers 
one considers closed loops on the lattice whose direction reverses at passive nodes h G Ad\A. Passive nodes (open 
circles (o) in Table |T] and Fig. are never crossed and reverse the direction of a path. Paths on the other hand do 
not reverse direction at active nodes. 

The four shortest loops on A^ of this kind are, 


Ca[n; := [n + A^, n, n - A^, n, n + Ap, n, n - A^, n, n + Ap] ; sgn[Ca] = e{pLiypa) , 

Cb[n; pvpa] := [n, n + A„, n + A„ + Ap, n - Ap, n, n - Ap, n + A„ + A^, n + A^, n] ; sgn[C{,] = £{pvpa) , 

Cc[n; pv] ■■= [n + Ap, n, n - A^, n, n + Ap] , 

Cd[n; pvpa] := [h, h + A„, h + A„ + Ap, n - Ap, n, h - Ap, h - Ap - A^,, h + A„, h] ; sgn[Cd] = £{pvpa) . (42) 


Using the assignments of Table. |Ij these closed loops correspond precisely to the densities of Eqs. (39) and ( |40| . 
Reflection at passive nodes is required by the local ^^{l) symmetry and ensures that spinors combine with complex 
conjugate spinors to anti-hermitian (null) lat-frames. Introducing the sign of a loop as in Eq.(42), the lattice action 
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of Eq.(411 for the duoverse can be succinctly written, 


*'S'hp°- = sgn[Ct]ViCb) - XI sgn[Ca]V(Ca) 


Cb 


(43) 


where the sums extend over the elementary loops of Eq.(42) on Expressions for the densities V{Ca[n', fJ-vpcr]) = 


and V{Cb[n] piypa]) = Vapui^jn + Aj) are given by Eq.(39a) and Eq.(|39b[) respectively. The imaginary 


parts of the expressions in Eq.(40bl and Eq.(40a) vanish. C^-loops with real coefficients give purely imaginary 


contributions to the lattice action. 

Field assignments and reflection rules thus dictate the form of the lattice action of the duoverse. Every configuration 
of the duoverse corresponds to two causally disjoint manifolds whose events are labeled by the hypercubic sub-lattices 
A and A'. The events on each of these manifolds are causally related to events of the same manifold, but not to 
events on the other. Parallel transport on spatial links between one sub-manifold and the other is possible, but the 
matter fields of each causal component do not interact. Matter and energy of one causally connected component of 
the duoverse thus influences the curvature of the other, but there apparently is no causal communication between the 
two manifolds. 


V. MAKINGS 010 CAUSAL MANIFOLDS 

Although causal manifolds are triangulated by null lat-frames as described in Sect0 not every lattice configuration 
corresponds to a triangulated manifold. For a causal simplicial complex the vertices of every B-tetrahedron must lie 
on the backward light cone of a future event. This condition is the same for uni- and duo-verses, since the two 
manifolds corresponding to A and A' lattices can be reconstructed independently. For definiteness we here consider 
the reconstruction of the A manifold. 


A. Consistency Conditions of the Simplicial Complex 

Necessary and sufficient conditions for constructing the triangulated manifold are obtained by considering the 4- 
simplexes of the pure complex. The complex is composed of two types of simplexes (which could be viewed as being 
local charts of an atlas), ch(n) and ch{n), 

ch(n) = {n, n ^ = 1,... 4} and c/i(n) = {n, n —/it;/r = 1,... 4} . (44) 

As schematically indicated in Fig. these simplexes are composed of an apex event and four causally related events 
on its forward or backward light cone respectively. We will refer to these special Minkowski simplexes with four 
light-like and six spatial edges as forward (backward) null-simplexej^ 

Four of the faces of such a null-simplex of Minkowski space are tetrahedrons with 3 spatial and 3 null edges, the 
edges of the fifth (F- or B-) tetrahedron are all spatial. Any two of these simplexes (charts) have at most two events 
(and thus at most one length) in common. The non-trivial intersections of such simplexes are (with p ^ v), 


ch(n) n ch(n + p) 

={n + p} 


ch{n) n c/i(n -|- p) 

={n,n-|-^}. 

with el^in) = el^{n + p) = 0 

ch{n) n c/i(n + p — v) 

={n-h/t} 


ch{n) n ch{n) 

={n} 


ch{n) n c/i(n -|- 2p) 



ch{n) n ch{n + p + v) 

={n -1- n -1- v}, 

with =Il^{n + p + u) 


Lyrics to the song ’’The Makings of You” by Curtis Mayfield: Add a little sugar, honeysuckle/ And a great big expression of happiness/ 
Boy, you couldn’t miss/ With a dozen roses/ Such will astound you/ The joy of children laughing around you/ These are the makings 
of you/ It is true, the makings of you/ The righteous way to go/ Little one would know/ Or believe if I told them so/ You’re second 
to none/ The love of all mankind/ Should reflect some sign of these words/ I’ve tried to recite/ They’re close but not quite/ Almost 
impossible to do/ Reciting the makings of you/ 

The corresponding 2-dimensional road atlas would be composed of charts that include just three cities, with one-way roads connecting a 
central city to the other two. The central city of any chart also is the central city of one other chart with reversed one-way roads (these 
two charts could be labeled past and future). In the universe every city is a central city on some chart. In the atlas of the duoverse 
peripheral cities (passive nodes) never are central ones on any chart. Every city is at the intersection of two one-way roads in both 
atlases. 
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FIG. 5. Color online. Schematic depiction of the elementary null-simplexes ch{n') and ch{n) forming the simplicial complex 
of the duoverse with passive (o) and active (•) nodes. In the universe all nodes are active. The four arrows of each simplex 
represent displacements to nodes on the forward (c/i(n')) and backward (c/i(n)) light cone of the active nodes n' and n. They 
are specified by the forward and backward lat-frames i5^(n') and Efi,{n) respectively. These illuminate F- and B-tetrahedrons 
whose 6 spatial lengths are shown in yellow and blue. Dashed lines connect nodes that are to be identified in the pure simplicial 
complex. The corresponding common edge of the F-and B-tetrahedron has spatial length I'i2(n — Ai — A2) = £i2(n). This is 
the geometrical content of Eq.(461. The four nodes (n, n — Ai, n — A2, n = n — Ai — A2 form a plaquette of the topologically 
hypercubic lattice with ^12(n — Ai — A2) = £i2{'n) on the diagonal. Note that the backward null lat-frames Efj, are redun dan t 
(and not included as variables of the model): if the 6 spatial lengths of the B-tetrahedron satisfy the inequalities of Eq.(491, 
the relative location of the apex n of the oriented backward null-simplex is uniquely determined. 


Simplexes (charts) of the duoverse include a single active node (n) and four passive ones.The first two intersections 
of Eq.(45) can be ignored in this case because they do not occur in the simplicial complex of the duoverse (i.e. are 
not in the atlas). 

Introducing the backward null lat-frame {£'^(n), /r = 1, ..., 4} on the backward light cone of an event n in the same 
manner as the forward null lat-frame {E^{n), /i = 1, ..., 4} was defined on the forward light cone of n, the simplexes 
(charts) can be consistently glued at their common overlap only if 


linin') := -2Ei^{n') ■ E^{n') = - 2 £;^(n) • Ej,(n) := , with n' = n - /r - 12 for all /r, 1 / and n . 


fLU \ 


(46) 


The same conclusion is reached by examining the typical /iJ 2 -plaquette of the A lattice shown in Fig. If the 
lattice conhguration is to represent a triangulated manifold, the events labeled hy n — fj, and n — v need to be on the 
backward light cone of n and on the forward light cone of the event labeled by n' = n — fa —v. Note that n —p, = xi! + v 
and n — 1 / = n' -b /r are the only nodes of the lattice in the intersection of these light cones. 

Including the case = i', Eq.(46l altogether imposes 10 independent S'L(2, C)-invariant constraints per site on 
the backward lat-frames. 

Equality of the time-like separation between events n and n' (see Fig. ) as viewed from inertial systems 
associated with n — fj, and n — v similarly requires that, 

0 = i?y(n — v) ■ E^{n — v) — E^{n — fj.) ■ E^{n — fj.) for all fx, v and n . (47) 


The se a re 6 additional independent linear constraints per site. The altogether 16 constraints per site of Eq.(46) and 
Eq.(47l offset the 16 degrees of freedom per site of the E’s. No additional dynamical degrees of freedom have been 
introduced and one thus should be able to construct a TLT that constrains lattice configurations to those representing 
triangulated causal manifolds. 

Eq.(46l as well as Eq.([4^ are invariant under local S'L(2, C) transformations of all the lat-frames, but Eq.(46) is 
invariant under independent SL{2, C) transformation of the backward and forward lat-frames at each active site. For 
a given configuration of forward null vectors {i?^(n); [n, /r] S A} the constraints (46) imply that. 


{0 < - p,-v) = -2E^{n - pL-v)- Ei,{n - p,-n), for 1 < /r < < 4} , 


(48) 


are the spatial lengths of the 6 edges of a B-tetrahedron whose four vertices are on the backward light-cone of the 
event n. 

The six additional constraints of Eq. (0 on the other hand ap par ently only serve to localize the additional 
S'L(2, C)-invariance E^{r\) —)■ g{n)Ef^{n)g''^(n) of the solution to Eq.(46). The ten S'L(2, C)-invariant constraints of 
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n 



n 



FIG. 6. Color online. Constraints of a triangulated manifold. Plaquettes of A in the duoverse are of type a), with two spatially 
separate passive (o) nodes and two time-like separate active (•) nodes or of type b), with two spatially separate active nodes 
and two time-like separate passive ones. Forward (backward) null lat-frames of the configuration are shown as red (green) 
arrows that point away from an active vertex and toward a passive one. In the universe all nodes are active. To construct a 
consistent simplicial complex, the spatial separation between the events labeled by the nodes n — fi and n — j/ must be the 
same for local inertial systems at n and n — ^ — u. The time-like separation between events represented by passive nodes 
n and n' also is the same for local inertial systems at the active nodes n — and n — v. In a) the center-event of a plaquette 
at the midpoint of the spatial diagonal is marked by a cross (x) and its time-like separation to the two active nodes indicated 
by solid (black) arrows (see Sect VIII. 


Eq. (46 1 in fact suffice to (uniquely) construct the whole complex. The geometry of this complex does not depend on 
the additional «S'L(2, C) degrees of freedom of the backward lat-frames. 

As in the GPS P5H?7] . the apex n can be reconstructed from the vertices of the B-tetrahedron of ch{n). The 
discussion of light-cones in Sect IV implies that this is possible only if the 6 spatial lengths ^^i/(n) = — ^ — v) > 0 


of Eq.(48l satisfy the inequalities, 


?L2(n)44(n) < TL3(n)£24(n) + luin^sin) , 

43(n)£24(n) < ?i2(n)£34(n) , 

£i4(n)£23(n) < TL 2 (n)£ 34 (n) -|- £ 13 ( 11 )^ 24 ( 11 ) , for all n . (49) 


Since all the lengths are positive, we again can combine the triangle inequalities of Eq.(49) to the single requirement 
that. 


det (n) < 0 , for all n . 


(50) 


The inequalities of Eq.(49) are 5^(2, C)-invariant and select configurations of null lat-frames that represent trian¬ 
gulated causal manifolds. Due to Eq.(50) the consistency conditions are the requirement that the 4-volume of every 


backward simplex V(n) := y —det£^^(n) be real. 

These quasi-local but non-linear consistency conditions on the lattice configuration not only are reasonable and 
necessary but also are sufficient to reconstruct the oriented backward lat-frame up to SL(2, C) transformation from 
the set of lengths {£^;,(n)} (see Appendix]^. As outlined in Sect|^ one then can reconstruct the simplicial complex 
of the triangulated manifold uniquely. The conditions of Eq.([50|) in some sense are ” integrability conditions” of the 
causal lattice. 


B. The Manifold TLT 


We here construct a local TLT whose partition function vanishes when Eq. (491 is violated and is a (non-vanishing) 


constant otherwise. The TLT partition function evidently is proportional to the product of Heaviside functions that 
enforce the inequalities of Eq.(49). We wish to have a local integral representation of it. Many such representations 
exist. In Appendix]^ an equivariant BRST construction is used to obtain a local TLT in terms of the lat-frames that 
ensures Eq.(|4^. Unfortunately this local TLT depends explicitly on the backward lat-frames and more than doubles 
the number of lattice variables. One also can simply enforce the single inequality of Eq.(50), which is of rather high 
order in the lattice variables. 
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Here we construct a simpler TLT based on the spinor formulation of the lattice model. It requires just two new 
variables (spinor phases) per site to enforce Eq. (491 with a local lattice action of the same scaling dimension as the 
cosmological term. 


Consider therefore the representation of backward null lat-frames E^^{n) by spinors 


(n) ■ 


* \ A.B' 


(51) 


The only difference to representation of forward lat-frames by spinors in Eq.(25) is the minus sign in Eq.(51). It 
implies that —iTrif^(n) = —(IC/Kn)]^ -b |5^(n)p) < 0. E^(r\) thus is on the backward light cone. 

Events labeled by n — /r and n — v are spatially separated by, 


and Eq.(46) demands that, 


V(n) = \din)eAB^!^in)\ 


^L(n) = ilM = 


(52) 


(53) 


A configuration of forward spinors {^^(n)} thus is compatible with a triangulated manifold only if real phases = 

(pvfiin) and backward spinors ^^(n) exist for which. 


(n) = := - n-v) for all . 


(54) 


The skew-symmetric matrix (n) defined by the right-hand-side of Eq.(541 is described by spinors only if its Pfafhan 
vanishej^ A configuration therefore represents a triangulated manifold only if a set of phases {(/J^y(n) = (p„f^{n) G R} 
exists for which 


0 = Pf(/(n)) = /i2(n)/34(n) -b /i3(n)/42(n) -b /i4(n)/23(n) for all n , 


(55) 


where /^iy(n) is given by Eq.(54|. Note that this formulation of the consistency condition assures the existence of 
backward spinors without explicitly constructing them. 

The phases can be found and Eq.(55) satisfied only if the magnitudes of the three complex numbers. 


a(n) 

6(n) 

c(n) 


= |/i 2 (n - Ai - A2)/34(n - A 3 - A4)| = £i2(n)£34(n) 

= |/i3(n - Ai - A3)/42(n - A 2 - A4)| = £i3(n)£24(n) 

= |/i4(n - Ai - A4)/23(n - A 2 - As)] = fL4(n)£23(n) 


(56) 


form the sides of a triangle. To construct the pure simplicial complex of a triangulated causal mani fold , the spatial 
lengths ^^i/(n) = |/^,,(n — A^ — Aj,)| = — ^v) thus must satisfy the inequalities of Eq.(|49|. 

Suppose that for a conhguration of forward spinors a set of phases can be found for which Eq.|^) is satisfied. 
We first demonstrate that a physically equivalent configuration of forward spinors exists in this case that satisfies 
Eq.(|55[) with vanishing phases. 


(57) 


The U^{1) invariance of physical observables (and lat-frames) of Eq.(28) implies the equivalence of phases, 

= (/ 7 ^^(n) - - fj. - 1^) - i)^{n - fj. - 1^) . 


For Pf[/(n)] = 0 the overall phase is irrelevant and we need only construct a representative set of spinors for which 
the phases of the three terms of the Pfafhan coincide, that is 


7512 (n) + 7 ’ 34 (n) = </ 5 i 3 (n) + 7 ^ 42 (n) = ■ 


(58) 


Given the initial set of phases for which Eq.(55l is satisfied, the 2N linear constraints of Eq.(58l on 

the AN phases ^;x(n) in general allow for an infinite set of solutions. To have a more manageable set we restrict 


14 j£ pf[j(n)] = 0, /^iy(n) = 7 (X^yi/ — XiyY^) where the two linearly independent vectors and span the two-dimensional kernel of 
the dual tensor ^ ^£(pupcF)fp(j{jY) and 7 is a complex number. A possible set of spinors then is ^/i(n) = Yp,)- 
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to V^i(n) = —' 02 (n) and ipai^) = —■*/'4('^)- Defining ^/’±(n) := V’i(n) ± i’si'n.) = —^’ 2 ( 11 ) =F V'4('^)i the conditions of 
Eq.(j^ then decouple into two equations, 

ip+in - A 2 - A 4 ) - V'+(n - Ai - A 3 ) = (/ 3 i 2 (n) + - <Pi 3 (n) - (^ 24 ( 11 ) mod 2t: 

^/;_(n--A 2 - A 3 ) - ?/)_(n - Ai - A 4 ) = (/ 3 i 2 (n) + (/ 334 (n) - (pi 4 (n) - (p 23 (n) mod 27r , (59) 


that relate the angles at the endpoints of two diagonals of a hypercube. With appropriate boundary conditions 
Eq.(59l uniquely determines the angles tj}±{n), mod 27r. The angles ^/’i(n) = (^/i+(n) +^_(n))/2 and '02(n) = 
('!/'+(n) — ijj-{n))/2 then are known modulo tt and the new spinors {^(j(n),/i = 1,... ,4} are determined up to sign. 


Let us therefore attempt to construct the TLT that enforces the manifold condition by imposing Pf[/(n)] = 0 at all 
nodes n as a C/^(l) gauge condition on the spinors. This ’’gauge” conditions has a solution only if the corresponding 
configuration of forward null lat-frames represents a triangulated manifolj^ 

We introduce ghosts c± and a corresponding BRST-doublet of complex anti-ghosts c, c* as well as complex Lagrange 
multiplier fields 77 , rj*. The nilpotent BRST-variation of the spinor phases ^p± and of these additional fields is, 


s'0±(n) = c±(n), sc±(n) = 0 sc(n) = 77(n), s? 7 (n) = 0 sc*(n) = r]*{n), sr]*{n) = 0 , (60) 

with a trivial extension to other fields. We therefore consider the TLT with the BRST-exact Lagrangian, 


.a, 


CTLT{n) = s[c*(n)Pf[/'(n)] -h c(n)Pf[/'*(n)] -h {n)c{n) + r]{n)c*{n))] 


= 77*(n)Pf[/'(n)] 77(n)Pf[/'(n)] iar]{n)ri* {n)+ 


(61) 


-I- i 


(c(n) c*(n) 


,/i3(n)/42(n) /{4(n)/;^2(n), 


c+(n - Ai - A 3 ) - c+(n - A 2 - A 4 ) 
c_(n - Ai - A 4 ) - c_(n - A 2 - A 3 ) 


where we have introduced a gauge-parameter a > 0. /^(n) in Eq.(61) depends on the spinorial phases and is the 
skew-symmetric matrix with components. 


/i2(n) = -/2i(n) := /i2(n - Ai - A2) 

/34(n) = -fisin) := /34(n - A3 - A4) 

/( 3 (n) = -/^i(n) := - Ai - A3) 

fM = -/ 42 (n) := e+*’^+("-^=-^-)/ 24 (n - A2 - A4) 

fUn) = := - Ai - A4) 

/ 23 (n) = -/^2(n) := e+*’^-("-^-^^)/ 23 (n - A2 - A3) . ( 62 ) 


The partition function of this TLT is. 


TT / dc+{n)dc-{n)dc{n)dc*(n) f d'0_|_(n)dV'-(n) f dr]{n)dri*{n)e 

n d JSixSi J\R^ 

cx TT / d7/'+(n)ci7/>_(n)P(n) . 

n Six Si 




(63) 


The last expression is obtained by integrating over ghost and Lagrange multiplier helds and 


P(n) := ilm/( 3 (n)/'*(n)/' 4 (n)/-(n) = I det [ 

V/l3(n)/42(n) /(4(n)/'2(n) 


(64) 


In the limit a —)■ 0+, only spinor configurations that satisfy Pf[/'(n)] ^ 
of Eq.(63). This is the condition we wish to enforce, but one can show 
the Euler characteristic of the 2N-dimensional torus T 2 N = {Si x Si)^. 
vanishes [38l Ho] for any a > 0 and spinor configuration {^}. This TLT thus fails to constrain to configurations that 
satisfy the consistency condition. 


^ 0 contribute significantly to the integral 
[55H55] tha1p^ [£] is proportional to 
As defined in Eq.(63l therefore 


The forward null lat-frames do not depend on the spinor phases and are 17^(1) invariants. 

Viewing Pf[/'(.)] as a smooth 2N-dimensional real vector field on T 2 N = (Si X Si)^ with coordinates (i/i+(.), i/i—f.)), this is a consequence 
of the Poincare-Hopf theorem l39l . 
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Fortunately this can be rectified. We already know that any skew-symmetric complex matrix with vanishing Pfaffian 
may be written in terms of spinors. Configurations that contribute significantly to the integral of Eq.(63) in the limit 
a —>■ O’*" thus can be written in terms of backward spinors as, 


/L(n) =L{n)eABL (n) . 


For Pf[/'(n)] = 0, P(n), given by Eq.(64), therefore is just the (signed) four dimensional volume element, 

P(n) = — |TrE 3 (n)eE^(n)eE 4 (n)eEf (n) = det[E(n)] , 


(65) 


( 66 ) 


of the backward null lat-frames given by Eq.(51). The lengths f^y(n) = |/(,j,(n)| themselves do not depend on spinor 
phases. Positive and negative contributions to in Eq.(63) thus arise from different signs of the 4-volumes 

E(n) of configurations with the same set of spatial lengths Non-singular solutions to Pf[/'(n)] = 0 in fact 

come in pairs with P(n) of opposite sign and this leads to = 0. 

Orientable manifolds on the other hand ought to be triangulated by spinor configurations for which V (n) > 0 at 
all sites n. It thus is tempting to modify the TLT to. 


lim n 

O'—>-0+ ^ 


a 


-1 


/ d'0+(n)fi^_(n)I/(n)0[I/(n)] e o|Pf[/(“)l| 
I Six Si 


(67) 


where 0[a;] is the Heaviside distribution. The question arises whether is a topological integral that does not 

depend on generic variations of the spinor configuration. This generally will not be true if the number of (pairs of) 
solutions depends on the configuration. 

We can assure ourselves that this is not the case by explicit evaluation of Eq.(67). Note that P(n), defined in 
Eq.(64|, depends on the following differences of the angles in Eq.(59) only, 


6»(n) = i( 4 '+(n- A2 - A4) - 4 '+(n- Ai - A3) -H 4 '_(n- A2 - A3) - 4 '_(n- Ai - A4)) 
(/j(n) = 4'4.(n - A2 - A4) - 4 '+(n - Ai - A3) - 4 '_(n - A2 - A3) -|- 'I'_.(n - Ai - A4) . 


( 68 ) 


The partition function of Eq. (671 thus factorizes in the form. 


^^^^[e]=n^[/i2(n-Ai-A2)/34(n-A3-A4),/i3(n-Ai-A3)/42(n-A2-A4),/i4(n-Ai-A4)/23(n-A2-A3)] 


with, 


^277 p27Z 


T[a, 6, c] := lim — 

a—>-0+ TTQ: 


d0 / dip Im(5c*e*‘^)0[Im(6c*e*‘^)] exp — 


+ t, ede+v^/2) + c 


(69) 


(70) 


For a —>■ 0^ the two-dimensional integral of Eq.(70) may be evaluated semiclassically. The exponent vanishes only 
at (absolute) extrema that satisfy the triangle inequa ities. 


|a| < |6| -I- |c| A |6| < |c| -I- |a| A |c| < |a| -I- |6| 


If Eq.(71) holds, these extrema correspond to angles 6,p that satisfy, 


= \b\'^ + |c|" -f 2Re bc*e^‘^ , = 


(71) 


(72) 


As Fig. illustrates, the exponent generically vanishes at pairs of absolute minimeP^ with angles p and p' related 
by be*= b*ce~^‘^ . Both extrema satisfy Eq.(72), but give the opposite sign for P(n) = Im bc*e''^. The Gaussian 
integral over fluctuations around these absolute minima furthermore are the same and proportional to a/|Im be*e^'^\ 
The total contribution to the partition function for a —> 0^ thus gives the vanishing Euler characteristic of the 2-torus 
if one sums over both minima. Taking only solutions with positive volume in Eq.(70l on the other hand leads to a 


If a is taken to be real, the two extrema correspond to letting {b,c) —^ (b*,c*). 
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FIG. 7. The pair of generic extrema with a + b + c = 0 for given magnitudes of the complex numbers b and c. The magnitudes 
|a|, |fo|, |c| satisfy the triangle inequalities of Eq.(|71[). Both solutions coincide only if |a| = ||fo| ± |c||. 


non-vanishing partition function that generically does not depend on a, b and c as long as these satisfy the inequalities 
of Eq.([7T|). T[a, b, c] on the other hand vanishes exponentially in the limit a —>■ 0 - 1 - when Eq.([7T|) is violated. T[a, 6 , c] 
thus is an integral representation of the distribution, 

T(a, b, c) = 0(|a| -h | 6 | - |c|) 0(|6| -f |c| - |a|) 0(|c| -h |a| - | 6 |) = 0 ( 20 ^ 6 ^ 2b‘^c^ + -a^ -b'^ - c^) , (73) 


Where 0(a:) is the Heaviside function. It is a constant if a, b, c satisfy triangle inequalities and vanishes otherwise. 

The TLT partition function of Eq.(67) thus constrains the space of configurations to those satisfying Eq.(49) that 
correspond to triangulated causal manifolds. Note that is a topological lattice integral in the limit a —)■ 0+ 

only. Eor a > 0, the integral does depend on the configuration. However , the limit a —> 0^ appears to be smooth 
and well-defined for |ap > a. Regularization of the model in Sect VHIC ensures that this bound is satisfied. 


VI. PARTIAL LOCALIZATION TO AN 5t/(2)-STRUCTURE GROUP 


The action of Eq.(|4I|) is real and the weight of a configuration therefore oscillatory. Although S'hp is un¬ 


bounded, it is conceivable that configurations with very large classical action do not contribute to the integral, as in 
ordinary Fresnel integrals. 

Even if this is the case, the Lorentzian model with the lattice action of Eq.(4I) is not well defined because 


its SL{2, C) structure group is not compact. The infinite volume of this group formally cancels in the expectation 
value of invariant physical observables, but it prevents one from defining a finite generating function even for a 
lattice with a finite number of sites. Contrary to ordinary lattice gauge theories with compact structure group, the 
SL{2, C) structure group of this model has to be localized to a compact subgroup for the generating function to be 
defined. The localization to an S'{7(2)-invariant lattice theory (where only the subgroup of local spatial rotations 
remains free) in fact is unique. The following partial gauge fixing of lattice configurations in this sense does not suffer 
of a Gribov [miig ambiguity. 

Consider the local 5'17(2)-invariant Morse function constructed from the spinors of a single active node n. 




(74) 


where {w^(n)} is a set of non-negative weights. This Morse function is proportional to an average over the positive 
time components = —i^Tig{n)E^{r\)g\n) of the S'L(2, C)-transformed null lat-frame. It is bounded below and 
invariant under the SU{2) subgroup of spatial rotations. Considered as a function of ^(n) € S'L(2,C) for a given 
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spinor configuration, thus is a function on the coset space SL{2, C) /SU{2) only. Decomposing g G SL{2, C) 

into a hermitian and a unitary part, both of unit determinant, 

g = uh, with h = h) = , and u G SU{2) , (75) 

critical points of are characterized by, 


0 = . 


(76) 


The 3x3 Hessian matrix H of this Morse function is strictly positive and proportional to the identity, 

H{n) = (n) = lWi[g{n)] . 


(77) 


det[i7(n)] vanishes only if ^^[^(n)] = 0, that is, when the null lat-frame is singulaip^ The non-trivial solution to 
Eq.(76l thus is unique (modulo spatial SU{2) rotations) and the Euler characteristic x(^(2, C) /SU{2)) = xi.Hz) = 1. 
Note that this partial gauge eliminates one of the Lorentz components of a forward null lat-frame as a dynamical 
variable. 

Several partial localizations of this kind are of special interest. They differ only in the choice of weights rc;i(n) > 0 
in Eq.(74). 


• The most practical partial localization uses just two linearly independent] null vectors at each node. Setting 
Wi(n) = W 2 (n) = 0 and ^ 3 ( 11 ) = W 4 (n) > 0, the gauge condition of Eq.(|76f implies that E 4 (n) = —E 3 (n) and 
T 3 (n) = r 4 (n) = |ii' 3 (n)|, i.e. the corresponding two events appear simultaneous in this local inertial system. A 
spinorial solution to this condition is, 


^i(n) =-el*(n), el(n)=^3i*(n) . 


(78) 


It is unique up to 17(1) phase transformations of the spinors. Note that Eq.(78) is compatible with the choice 
made in Sectj^to solve the consistency constraints with il’sin) = Due to its simplicity, the partial 

localization of Eq. (|7^ will be used in much of the following. 


• It may be desirable to choose equal weights for three of the four components, i.e. rc 4 (n) = 0 and Wi(n) = 
W 2 (n) = rc 3 (n) > 0. The spatial components of three null vectors here form a triangle whose side-lengths are 
the temporal components of these null vectors. However, an explicit solution to Eq.(76) in terms of the spinors 
is quite involved. 

• Sometimes a more equitable gauge condition is preferable that weights all components equally wi(n) = W 2 (n) = 

r(; 3 (n) = u> 4 (n) . Eq.(76) then implies that the spatial parts of the four null-vectors of the lat-frame sum to 

zero. This gauge condition is used to reconstruct the backward lat-frame in Appendix [C| 


The weights w^(n) in Eq.(74) can be any positive S'T(2, C)-invariants of the fields. A judicious choice of these 


invariants may sometimes be advantageous (as’t Hooft gauges are in spontaneously broken gauge theories). In 
the next section we have occasion to consider path-dependent gauges where the weights u)^(n) depend on n in 
a manner that simplifies the computation of a path-dependent quantity, such as the proper time. 

The partial localization of 5'L(2, C)to the SU(2) subgroup at each site of the lattice may also be achieved by an 
equivariant BRST construction [35113S] • Since the gauge condition here is ultra-local, the construction here is rather 
trivial, but it illustrates the general principle. The nilpotent BRST variation of the spinors, ghosts (c and uj) and 
auxiliary fields {h and c) is. 


= ^(c +w) • , 

sc = UJ X c , 
sc = b + uj X c , 


SUJ = ^C X C+ X OJ 
sb = uj X b + ^c X (c X c^ , 


(79) 


The invariant regularization of the theory in 


Sect 


Vine 


implies that det B(n) 


det (_E“(n)) > 0 for all n. 
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Note that the w-ghosts here generate SU{2) transformations. The partial localization then is implemented by extend¬ 
ing the S'L(2, C)-invariant lattice action by the 5'C/(2)-invariant and BRST-exact part, 

Sgf = s^c(n) • I y (6(n) -h/cc(n) x c(n)) + ^ w^(n)C^(n)CT^^(n) j 

n \ fi / 

= (*7«:c(n) xc(n)-t^W;,(n)^^(n)CT^^(n) I -f c(n) • c(n) ^ n;^(n)^^(n)ie^(n)-f 


+ y Kn) • b{n) - y (c(n) x c(n)) • (c(n) x c(n)) 


'‘Yl [2 (^n) + '«c(n) X c(n) - ^ ^ n;^(n)^^(n)CTC^(n) j -t- ^ | ^ (n)^y („) j 


■ in (c(n) X c(n)) • ^ zi;^(n)e;;(n)aeM(n) - *c(n) • c(n) ^ w^in)^l{n)l^^{n) + 


+ Y ^ \ c(n) X ?(n)) • (c(n) x c(n)) 


(80) 


where n and 7 > 0 are gauge parameters. Observables and the extended lattice action, do not depend on the w ghosts 
that generate local SU(2) variations. The Grassmann integrals of w therefore can be saturated and the equivariant 
BRST variation is obtained by formally setting w = 0 in Eq.(79l, 


SeC/J 2 ’ 

SeC= 0 , 

SrC = b , 


Sed = 


6 = ic X (c X c) . 


(81) 


The equivariant BRST variation Sg is nilpotent on SU(2)-invariant functionals only: form ally Sg generates SU(2) 


variations with gauge parameters ^c(n) x c(n). Note that the quartic ghost interaction of Eq.(80) vanishes for k = ±1 


or 7 = 0 only. Since Sgf is ultra-local, this part of the action can be absorbed in the integration measure at each site. 
The fact that the expectation of S'L(2, C)-invariant observables does not depend on the gauge parameters implies 
that the limit 7 —> 0^ is smooth. One recovers the gauge condition of Eq.(76) in this limit and the Grassmann 


integral over the ghost fields provides the factor det[i7(n)Jin the local bosonic measure of the spinors, where the 
Hessian H{n) is given by Eq. ^GF shows that Eq.([7^ is just one of many ways to localize SL{2,C) to the 
compact SU{2) subgroup. 


VII. PROPER TIMES OF CAUSAL PATHS 

A basic quantity of interest is the geodesic distance between two nodes of this lattice. Gontrary to ordinary lattice 
gauge theory, this distance depends on the lattice configuration and is not specified a priori. The proper time of 
a causal (time-like) path between two nodes depends on the configuration and correlation functions will have to be 
conditional on this proper time. 

In Sect. 1^ we saw that lattice configurations which satisfy the inequalities of Eq.([4^ can be interpreted as 
triangulated causal manifolds (or of two disjoint causal manifolds in the case of the duoverse). Here we define 
geodesic distances on these manifolds and in particular define the proper time of a causal path between two nodes for 
any consistent lattice configuration. 

Two nodes n and n' of the A- (or of the A'-) lattice are (uniquely) related by, 

n' = n -I- ^ . (82) 

The four integers {h^} give the following causal relations between the nodes n and n', 

hfj. > 0, for all /r = 1 ,... , 4 , n' > n; n' is in the future of n 

h)j. S: 0, for all /r = 1 , ..., 4 , n' < n; n' is in the past of n 

otherwise , n' / n; n' and n are not causally related . (83) 
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A lattice path of depth L from node ng to node is a list of contiguous nodes, 


‘Pl = [no,ni,...,ni] = {no,nj = ± for i = 1,..., L} . 


(84) 


Since the separation between any two adjacent events of the null lattice is light-like, the proper time of a path like 
vanishes. To allow for paths that are time-like, additional events of the manifold must be considered. It is natural 
to include events that bisect the spatial diagonals of plaquettes. In Fig. these are marked by a cross (x). These 
additional (centered) events of a plaquette based at the active node n will be denoted by n. Note that the plaquette 
is uniquely determined by the events on the path before and after n. The temporal and spatial separation of these 
center-events to other events of the same plaquette are readily computed in terms of the null lat-frames of a single 
inertial system. 

Consider a /ii^-plaquette based at the active node n with ’’center” ri. Viewed from an inertial system with origi n at 
the active node n, the event proper time to the event n is As = |A^(n) • Ei,{n) = i£^^(n) > 0. Due to Eq.(46) 

the proper time between ri and if = n + fi + v on the same plaquette. determined in the inertial system at n', is the 
same As' = iA^(n) • i?y(n) = = As . Note that the proper time 2As = between n and n' is largest 

for a path that passes through the central node n of the diagonal. The spatial separation of the event at n from the 
events labeled by n + /r and n + v alscplis \£fi^(ii). 


Extending Eq.(84| to include paths ‘Pl of depth L between two (active) nodes n and n' that can pass through 
central nodes, we consider contiguous lists of events of the form. 


= {no,...,nj_i,rij, nj+i,...,nfc_i,rife,nfc+i,...,nL} . 


(85) 


Note that any center-event rife in this list is flanked by two (active) nodes of the same plaquette. We can select the 
causal paths from an (active) node ng to an active node > rig by demanding that each of its increments either be 
light- or time-like in the forward direction. A causal path thus is a list of L contiguous events of the forirp^ 


Cl = {ng,... ,nj_i, n^, Hj+i,..., nfc_i, rifc, n^+i,..., iil} , 


( 86 ) 


with time- or light-like positive increments only. Every point of a causal path thus is in the immediate future of the 
previous one. The adjacent active events determine the plaquette {fJ-r') on which a center-event rife is located. The 
proper time s(Cl) of a causal path is the sum of the proper times of its increments. 


s{€l) ■= ^ 2^-\Ef,{ii.k-i) ■ E^{ii.k-i) = ^ |^i.(nfc_i)| > 0 , with rifc+i = nfc_i -f -I- , 


SfcGCi, 


(87) 




The sum is over center-nodes of the path only, since all other increments are null and do not contribute to the proper 
time. Note that the proper time defined by Eq.(87) is an observable composed of basic invariants of the lattice. 


The consistency condition of Eq. (46) has to be satisfied, but the backward lat-frames of the conhguration need not 


be reconstructed to define this proper time. That Eq.(87) is a physically sensible definition of the proper time of a 


causal path can be seen by choosing inertial systems (a gauge), in which the spatially separate events of any traversed 
plaquette are simultaneous (we proved in Sect VI that this path-dependent choice of gauge is possible and unique up 


to spatial rotations). In the corresponding inertial systems, the proper time between the two active nodes also gives 
the spatial extent of the plaquette. Eq.(87) merely defines the proper time of the whole path as the sum of 

these incrementj^ Fig. [^depicts some causal paths of depth L = 4 and their proper times. 

For any particular consistent causal lattice configuration one may select the lattice path with the longest proper 
time and define the geodesic distance between two causally related nodes as. 


d(n, n') := max s(£l) , 

:n-An' 


( 88 ) 


where L = /^^ is the depth of all causal paths that relate n to n^ The number of causal paths that connect a 

given pair of (active) nodes is finite but can be rather large: the number of causal paths with vanishing proper time 
already is L\/{hi\h 2 \hs\h 4 \). In general d(n,n') = 0 only if n' = n + but the geodesic distance between some 


Although the 51/(2, C) invariant square of the distances have opposite sign, proper times and spatial lengths here are both defined 
positive. 

Every causal path between two active nodes n and > n has the same depth L = hf_i . 

Note in this context that for a 34-plaquette in an inertial system at n that satisfies Eq.( |78[ l, |/ 34 (n)| = = ^3 = T 4 . 
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FIG. 8. (Color online) Nodes of the duoverse related to n by causal paths in the /ij/-plane. Several causal paths of depth L = 4 
through active nodes (•), passive nodes (o) and central nodes (x) are shown in different colors. Solid sections of a path do 
not contribute to its proper time. Dashed sections of a path through a central node contribute to the proper time, where 
(■iiv > 0 is the spatial separation between the two passive nodes of the same plaquette. In the universe all nodes are active and 
there is a central node on each plaquette. 


events may be quite small (of Planck size) even for a large number of steps L. From a macroscopic point of view such 
events essentially are on light-like geodesics. 

These considerations apply to a given causal configuration only. Since the weight of the Lorentzian lattice is 

complex, the weighted proper time of a given path between two causally related nodes in general could be complex 
even though it is a positive real quantity for any particular configuration. Stationary phase approximation suggests 
that the average proper time is (almost) real for sufficiently deep paths, but could be complex for paths of a few 
steps. The geodesic distance between close-by nodes depends on the configuration and a physical interpretation of 
correlation functions between two local observables (say P(n) and P(n')) such as 

(P(n)l/(n')) , (89) 

is hardly possible if the geodesic distance between the lattice points n and n' fluctuates. One perhaps instead should 
consider conditional amplitudes like. 


(y(n)F(n'))d(„^„')=T 


(90) 


Here the correlator is the (complex) weighted average with respect to configurations whose geodesic distance between 
the two lattice vertices n and n' is given by r. One then could ask how this amplitude depends on r > 0. The 
conditional amplitude of Eq.(90) thus gives the correlation between the local 4-volumes of a subset of all causal 
configurations. Since the theory is not invariant under translations, the conditioning of amplitudes in general will 
have to be more severe. It may in particular be necessary to fix the geodesic distance of an event to the Big Bang, 
because correlators like that in Eq.(90) almost surely depend on the epoch. 


Eor nodes n and n' related by causal paths of great depth L, the conditioning in Eq.(90) may not be that important. 


because the geodesic distance in this case is relatively well defined and should not differ much from the leading semi- 
classical estimate. Conditional and unconditional correlators coincide in this approximation. One thus perhaps can 
compute the asymptotic large-distance behavior of amplitudes in the conventional fashion. However, only conditioned 
correlators may be well defined at short distances on a finite lattice. 


VIII. INTEGRATION MEASURE, REGULARIZATION AND ORIENTATION 

To compute expectation values of SL{2, C) -invariant observables, an integration measure for the spinors and trans¬ 
port matrices must be specified. Contrary to ordinary lattice gauge theories with fixed lattice spacing, the short- 
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distance behavior of the present lattice model in addition needs to be regularized in a coordinate invariant geometrical 
fashion. The proper measure satisfies both objectives and in addition ensures that configurations of this lattice model 
correspond to oriented triangulated manifolds. 


A. Polar Parametrization of Spinors 

In a ’’spherical” parametrization of the spinors by their magnitude r and (real) angles 0 < ijj^ < 27r, 0 < < 2tt, 

0 < < TT, 




/ e cos(ej,(n)/2) \ 

sin{e^{n)/2) J ^ 


T^(n) is the only non-compact variable, 


(91) 


0 < r^(n) = ^ |^;^(n )|2 = -iTiEf^in) . 


A^l 


(92) 


f7^(l)-invariance of physical observables implies that they do not depend on the phase ^/'/^(n) of the spinors. In 
spherical coordinates the singular anti-hermitian matrix E^{n) thus is parameterized by its temporal extent, r^, and 
the three-dimensional unit vector, 


V(n) = (sin 6 l^(n)cos(p^(n),sin 6 »^(n)sm(/j^(n),cos 6 i^(n)) , 

which specifies the direction of the null-ray in a local inertial system at n. 

In the parametrization of Eq.(911, the forward null lat-frames may be written. 


= f^^(n) (g) (n) = iT^(n)P^(n) , 


where the projector is the matrix. 


(93) 


(94) 


P^(n) = P( 7 >(n)) :=-(!- P 


,(n)) = i( 


l-|-cos(e,j(n)) e sin(e,i(n)) 

sin(6l,j(n)) l-cos(6(,i(n)) 


)^ 


(95) 


Using Eqs. (94) and (95) in Eq.(39a), the 4-volume element U(n) of the 5-simplex formed by the four null vectors 
of a node is given by, 


1 


ri 1 


U(n) = —TiT 2 r 3 T 4 det ( ^3 M = Ti(n)r 2 (n)T 3 (n)r 4 (n)y(n) 

V U 1 / n 


(96) 


with 


T’i-r4 

r2-r4 

r3-f4 


U(n) = — det ( ? 2 -U ) = X! £(MJ^/Ocr)TrP^(n)eP,.(n)ePp(n)eP,^(n)e . 


I 

Is 


(97) 


flipper 


The 4-volume thus is proportional to the three-dimensional volume of a tetrahedron with vertices at the unit vectors 
U)^ 2)^3 and fq. It vanishes only when these four points are coplanar and changes sign when a vertex is reflected 
through the plane of the other three. 

In a system where Eq.(78) holds, the dependence on T 4 and fq can be eliminated. The 4-volume in this case 
assumes the form, 


VGF(n) 


fi 1 

r2 1 


- = -n 

ra 1 I o i 


('niTTln')r.?fi 


(98) 


and is proportional to the three-volume spanned by the spatial components of Ei{n), E 2 {n) and £' 3 ( 11 ). 
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B. Localized Lattice Integration Measure 


Since transport matrices and lat-frames are invariants, the diffeomorphism group does not constrain the integration 
measure. However, invariance under the 5'L(2, C) structure group and the local C/^(l) to a large extent dictates the 
local integration measure of this lattice model. The invariant integration measure for the SL{2, C) transport matrices 
is the Haar measure of SL{2, C). For matrices in the fundamental representation of SL{2, C) of the form, 

^”( 75 )’ (a, /3, 7 , (5) S and aS — = 1 , (99) 

the left and right invariant Haar measure is proportional to [43]. 

dfj,[U] = + /?7 — a6)dada*d/3dl3*djdj*dSdS* = dada*dl3dl3*d'jd'y*/\a\'^ . (100) 


For the purpose of analytic continuation of SL{2, C) to a compact SU(2)]^ x SU{2)fi it may be more convenient to 
parameterize U G SL{2, C) in terms of 5 compact Euler angles 0 < tp < 27r, 0 < 9,9 < tt,0 < < 4tt and a single 

non-compact variable 0 < pj < 00 , 


U — gO-ad/S ^cj.2612 gCT30/2 gCr 2 e /2 ^rJ3(t)/2 . jj^ _ ^-(730/2 ^-0^6/2 /2/2 ^-cji,612 g-o-3-0/2 ^ ( 101 ) 

The corresponding S'L(2, C)-invariant Haar measure is proportional to, 

d^\U\ = dip sin 9 dddcpdp) sin 9 dBdcp , (102) 

over the entire parameter domain 

Invariance under SL{2, C) transformations of Eq. 0 furthermore greatly restricts the local integration measure for 
the null-vectors E,^{n) and for the corresponding spinors. For the spherical representation of Eq.(94| the SL{2,C)- 
invariant measure at each lattice site is proportional to. 


d^[E,j] oc d'^E^, S~^{Ei, ■ E,^) oc T^^d'^r^, dfl{r,j) . 


(103) 


dH(fv) = sin{9,j)d9,jd(pi, here is the S'0(3)-invariant measure on S 2 in spherical coordinates and = 0[Li]dT^. 
The corresponding S'L(2, C)-invariant measure for the spinors includes the integration measure for the phase '0^. 
Local 17(1) invariance requires that. 




,d'^T,jdipydVt{r,y) 


where the domain of the three angles is 0 < ip^, < 27r, 0 < 9,^ < n and 0 < (fu < Stt. In the gauge of Eq.( 78) 
integration measure for the spinors at each node thus is of the form. 


(104) 
the local 




dtpi p{V) . 


(105) 


The local {7^(l)-invariance of observables was here exploited to set tp2 = — ipi and ipi = —tps and the d eterminant of 
the Hessian of Eq.(77) leads to the weight proportional to t| in Eq.(105). The measure of Eq.(105) distinguishes 
/i = 3 from p = 1,2, since the partial SL{2, C) gauge-fixing of Eq.(78) does not treat the lat-frames equitably. 

Requiring SL(2,C)- invariance thus determines the local integration measure up to an invariant density. p(V) > 0 
in Eq.( 105) is a positive semi-definite function of local SL{2, C) -invariants. All these local invariants are constructed 
from /^7In) and they include the 4-volume V (n) of Eq.([^. Note that p(V) in the presnt model generally will depend 
on the local 4-volume: in the first order formulation of Appendixj^ auxiliary 0-forms are introduced to linearize matter 
actions. Integrating over such non-dynamical fields leads to a factor of 1/ \JV (n) in the lattice measure at each site for 
each auxiliary degree of freedom. p{V) must compensate for the introduction of such non-dynamical and unphysical 
degrees of freedom. Note that the dependence of the classical lattice action on the cosmological constant could also 
be absorbed in p{V). 

Assuming that p{y < 0) = 0 resolves the sign ambiguity of the volume element and ensures that all lattice 
configurations correspond to oriented complexes. The next section exploits the fact that the critical limit of the 
model depends on the asymptotic behavior of p{V) for F —>■ 0 only. 


Formally the analytic continuation of SL(2,C)to compact SU{2)i^ X SU{2)ii corresponds to a change of the integration contour for ■ip 
from the positive real axis to the imaginary one, ip —> ±iip. However, this ”Wick rotation” of the high-dimensional integral must be 
carried out with care|44| for the contribution from the contour at infinity to vanish. Although important for numerical simulations, this 
article does not further pursue the analytic continuation of the lattice integrals. 
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The partially localized integration measure of Eq.(105) is manifestly invariant under the residual SU(2) structure 
group. The local rotational invariance of observables therefore can be exploited to completely localize the residual 
SU{2) structure group. This localization of a compact group is not compatible [5511101 with a BRST-construction, but 
one evidently can choose lat-frames in which the /r = 3 ray has spherical coordinate 83 = 0, with ip^ undetermined, 
and the fi = 2 ray defines the a:0-plane with 1^2 = 0. Renaming tpi —>■ ip, the local 4-volume of Eq.(98l in this 
complete gauge is given by, 


V'GF(n) = -Ti(n)r2(n)T|(n)sin6»i(n)sin6»2(n)sin(^(n) . 


(106) 


The local volume element V in this case is positive for 0 < (/? < tt and negative for tt < </? < 27r (since 0 < 81,82 < tt). 
One thus can restrict to oriented lattice configurations by the integration domain of p ! For oriented configurations, 
the integrals of the remaining spinor variables at each active node n of the fully localized model are, 





/ sin 02^6*2 / sin^id^i 
/o Jo 



where the spinors associated with a node are parameterized as. 


di!)! p(V'gf) 


?4 = v^e-^M0) 


= ; 6 = V^e 


— iip 


1 / cos (62/2) \ 
sin (62/2) ) 


6 =(:; 


cos(ei/2) ' 
sin(ei/2). 


(107) 


(108) 


Ignorinthe remaining spinor angles i/'i and ^3, 
5'L(2, C)-invariant degrees of freedom. 


Eq.(107) is an integral over the correct number (6 per site) of 


C. Invariant Regularization 


Although discrete, this lattice model is not necessarily regular at short distances and may not possess a critical 
limit. Configurations for which the local lattice volume E(n) is arbitrary small are not suppressed by triangulating 
the manifold and may dominate the lattice integral. A minimal coarseness must be imposed to avoid the associated 
UV-instability and obtain a UV-finite lattice model. The UV stability of the model thus is related to the behavior of 
the density p{V) at small V. Note that the local 4-volume element V (n) is invariant under reparametrization because 
the lat-frames are (see Eq. @). 

Since events are defined by the intersection of light-cones, one cannot demand that the local 4-volume be of fixed 
size. A viable model requires only that p(V) suppress configurations with small local 4-volumes E(n). In the critical 
limit when the number of lattice sites JV becomes large, while the total 4-volume of the universe remains constant, 
only the asymptotic behavior of p{V) for small V matters. I will assume that it is power-like 

p(E - 0) cx 0[R] E^(l-f C>(y)) , (109) 


with an a priori unknown global exponent 7 that could depend on the number of lattice sites N and on the cosmological 
constant A (and possibly other lattice couplings). 

In the complete gauge of Eq. (78), the asymptotic UV behavior of the density in Eq. (1091 leads to lattice integrals 
at each (active) node of the form, 

pOO nTZ pOO pOO pTT p7T p‘2'K p 

/ / {sinpydp / T 2 ''’^dT 2 / dri / {sin82)^^"'882 / {811181)^'^'^ddi / d-ips 

Jo Jo Jo Jo Jo Jo Jo Jo 


2 -k 


dipi , ( 110 ) 


where the parametrization of the spinors is given in Eq.(108|. 

For a lattice with N sites, the total 4-volume of the uni- (or duo-)verse will be constant on a curve 7(A; N). The 
intersection of these curves in the limit V —> 00 then formally yields the critical points 7*, A* of the continuum theory 
(if a critical limit exists). 


The integrals over tpi and 1/13 in 
not depend on these angles. 


Eq.l |l07| are absorbed by the TLT of Eq.( |67| l that enforces the consistency condition. Observables do 
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IX. THE STRONG COUPLING LIMIT A ^ oo 


It is interesting to consider the (naive) strong coupling limit of this model in which curvature terms of the action are 
simply neglected. The lattice action in this limit is ultra-local and the model has a number of symmetriej^ These 
enable us to define a finite lattice model and establish that 7 * = —2 in this strong coupling limit. At this critical 
value of the exponent 7 the 4-volume of the universe collapses for any lattice with a finite number of sites. 

Using Eq.(108) and the definition of Eq.(32a), the lattice integrals of Eq.(llO) at each active node n may be 


written in terms of the invariant lengths of Eq.(33) as, 


n 


Mt 


V^-^T{aAc) , 


with T defined in Eq.(|73[) and, 

a(n) = £i2(n)£34(n), 6(n) = ^i3(n)£24(n), c(n) = ■ 


Expressed by the lengths, the local invariant volume element of Eq.(39a) is. 


V (n) = det(£2^(n)) = {2a?}p' + 26^-I- 2c^a? — — b‘^ — c'^) 


4\1/2 


( 111 ) 

( 112 ) 

(113) 


Note that the triangle inequalities imposed by the T-distribution in Eq.(lll) arise automatically in the change of 
integration variables. They ensure that the volume in Eq.( |113 ) is real. The fact that the localized measure of 
Eq.(llO) can be expressed in terms of 5'L(2, C)-invariants indirectly confirms that the localization was unique. 

Since the cosmological part of the lattice action Ssc ~ ^(’^) ultra local, the generating function of the naive 

strong coupling limit decomposes into a product of independent generating functions for each site when consistency 
constraints are ignored. We restrict the configuration space to consistent complexes (that represent triangulated 
causal manifolds) by imposing the inequalities of Eq.(49), respectively Eq.(50). The generating function Zsc[^\l] 
in the naive strong coupling limit we are considering formally thus is. 


Zsc[^\l\ = 


nGA 


n 

\tl<V ‘ 


Mi 


yi 


-l-AiXV 


T(a,5,c)T(h, &,b) 


(114) 


where a, 6 and c are_given by the ^^i,(n) = .^^y(n — ^ — z/) as in Eq.(56). The T-distributions of Eq.(731 ensure that 
the integral in Eq.(114) is over configurations {(.^^(n)} that correspond to triangulated causal manifolds only. Note 
that the triangle inequalities of Eq. (731 actually depend on the squares of the lengths only, as might be expected 

since the continuum action in this limit depends on the metric gik only. _ 

For a hypercubic lattice with active sites and periodic boundary conditions, the variables of Eq.(112) and 
Eq.(56l are invariant under 3iV^ continuous scaling symmetries. These form an Abelian dilation group generated by, 


S^Hi 4 {n') = U 4 (n') for all n' = n-|-with g G Z and {i,j,k} = {1,2,3} 


where := A, 


Afc — Ai — A4 


(115) 


This definition of the generators may appear to single out the 4-th lattice direction, but in fact implies only that 
scales oppositely to ipa, with {fjL,v,p,a} = {1,2,3,4}. By definition if n and n' differ by a multiple 

of A^®). There thus are only 3A® independent generators on a periodic N'^ lattice and the number of these scaling 
symmetries is proportional to the number of sites on the 3-dimensional surface of a hyper-cubic lattice. If one ignores 
the consistency constraints, the number of independent scaling generators increases to 3A^ and is proportional to the 
number of lattice sites. _ 

'u Je 


One readily verihes that the integration measure 


n,/i.<L 


of Eq.(114) is invariant under the dilation group 


generated by Eq.(115). The periodic lattice thus is invariant under 3A® non-compact scaling symmetries and the 


generating function of Eq.(114) diverges. 


However, the expectation of observables that are invariant under these scaling symmetries, such as any lattice 4- 
volume, can be computed by localizing the integrals with respect to the scaling group. More interestingly, this dilation 


Although the evolution of our universe is believed to be dominated by the cosmological term long after the Big Bang, we here consider 
rather small lattices that hardly provide a faithful representation of the universe in its later stages. 
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FIG. 9. Conic lattice: an impressionist rendering of the conic part of the 4-dimensional universe that is causally connected 
to the Big Bang (the event labeled by (0,0,0,0)). [This figure does not represent results of a numerical simulation.] The four 
’’spinal” null rays on the light cone emanating from the Big-Bang event are depicted together with events on them. Spatial 
hyperplanes labeled by their depth L are indicated schematically. Note that these hyperplanes do not correspond to cosmic 
time because the vertices of the spatial tetrahedrons of a null-simplex generally are not simultaneous in any inertial system. 
The sketch presumes that these hyperplanes rapidly approach Ss-hyperspheres with increasing lattice depth. The connected 
part of the duoverse sketched here is rather small with a maximal depth L = 7, that is all its events can be reached in < 7 
light-like steps from the Big Bang. 


group generally is broken by boundary conditions at the ’’surface” of the lattice. For physical reasons, we would like to 
model the causally connected part of the universe by a conic four-dimensional lattice like the one sketched in Fig. 
This conic section is obtained by standing a hypercubic lattice on one of its corners and cutting it off at the spatial 
hyperplane with depth Lniaxl the apex at n = (0, 0,0,0) of the cone in this simple model represents the Big Bang 
(BB), all other events on the lattice being causally related to it. 

T he q uantities a,b,c of Eq.(112) are invariant under the scaling symmetries by construction and the 4-volume of 
Eq.(113) is as well. E^(n) furthermore is a quadratic form in a^,b^,c^ which in addition is invariant under a local 
SL(2,R) symmetry. Although this local SL(2,R) may be an accidental symmetry of the action and explicitly broken 
by curvature contributions, it also is a non-compact symmetry of the local integration measure without consistency 
constraints. If consistency constraints are imposed, only a subgroup of local SL(2,R) transformations at boundary 
nodes is preserved. _ 

In Appendix the scaling symmetries of Eq.(115) and the local SL(2,R) are examined in detail. It is found 
that these non-compact residual symmetries can be localized by imposing boundary conditions on the cone. Viewing 
these boundary conditions as a localization by the symmetry, one can adjust the integration measure so that the 
expectation of invariant observables does not depend on them. The boundary conditions of Eq.(EI7) effectively 
identify parts of faces of the global spatial tetrahedron of given depth L. Note that this identification turns the 
triangulated 3-dimensional spatial hyperplane of depth L of the hypercubic lattice we started with into a triangulated 
closed 3-manifold. _ 

Ignoring restriction of the integration domain imposed by T(a,b,c), the generating function Zgc of Eq.(114) 
factorizes in generating functions at each site. The local SL(2,R) and scaling symmetries in this case (see Appendix |E) 
allow one to deduce that. 


^sc|A;.|=(Vf) 


N(y+2) 


(116) 


where N is the number of sites of the lattice and z(j) for 7 > —2 is a finite integral that does not depend on A. Erom 
Eq.(116), the expected 4-volume of the lattice universe at strong coupling thus is given by, 


= - j A l..i:sc|A;7l = . 

when consistency constraints are ignored. The variance of the total 4-volume in this case is, 


>" Zso|A;7| = bT|V 


(iRUniv.^ 


SC 


4A 


- > 0 . 


(117) 


(118) 


The standard interpretation of Zgc as generating (here depth-) ordered vacu um e xpectation values Since the 
expectation of the variance of an hermitian operator ought to be positive, Eq.(118) implies that the 4-volume 
of Lorentzian space-time should be associated with an anti-hermitian operator V on the Hilbert space and that 
one measures eigenvalues of This is in keeping with anti-unitary time conjugation in ordinary space-time: 
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although V (n) changes sign, iV (n) is invariant under time reversal. This assignment is compatible with the fact that 
the generating function of the strong coupling limit turns into a real probability measure upon analytic continuation 
4iA ^ /3 > 0. 
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FIG. 10. (Color online) Numerical results for a consistent conic null lattice at strong coupling. The cosmological constant in 
this numerical study is purely imaginary /3 = 4iA = 1. 7 > —2 is the exponent of the density piV) oc V~' . a) Running averages 
of the invariant subvolumes {V(L)) up to depth L < Lmax = 7 of the conic lattice. {V(Univ.)) = (V(Lmax)) is the average 
total 4-volume of the conic null lattice, b) Distribution of the total 4-volume of the conic null lattice with Lmax = 7 for /3 = 1 
and 7 = —1.95. c) Dependence of the average 4-volume of the null lattice on the depth to the apex of the cone for /3 = 1 and 
7 = —1.95. Note that the 4-volume increases much slower than the number of nodes, Anodes oc L'^. d) Total 4-volume of the 
cone up to depth Lmax for /? = 1 but different exponents 7. The numerical results for 7 = —1.5, —1.8, —1.95 and a constrained 
linear fit to them are shown. The extrapolation to vanishing 4-volume at 7 = 7* ~ —2 is within the numerical accuracy of the 
simulations. 


Inclusion of the consistency conditions requires numerical simulation even in the naive strong coupling limit. The 
numerical results for a conic null-lattice with a maximal depth Lmax = 7 (corresponding to 330 nodes) for a purely 
imaginary cosmological constant f3 = Ai\ = 1 are presented in Fig. |10| For these simulations the residual non-compact 
’’surface” symmetries were all localized in the way described in Appendix]^ For ImA < 0 , the corresponding generating 
function is finite and a Metropolis-Hastings algorithm was used to obtain the shown results for /3 = 1 . Fig.[Io|l) shows 
that the 4-volume of the cone vanishes for 7 ^ —2. The constraints imposed by the consistency conditions therefore 
do not change the critical exponent 7 *. This is not too surprising, since the consistency constraints are irrelevant for 
the UV-behavior of the model. To see reasonably fast convergence on the finite lattice, we simulated the model for 
7 > —1.95. Fig. 10 r) shows expected volumes {V{L)) of the cone up to depth L, (V(7)) = ) corresponding 


to the volume of the null lattice cone itself. These are some of the invariants of this model. They all converge even 
for 7 = —1.95 ~ 7 * with the largest fluctuations in 14(0), the 4-volume of the simplex whose apex is the Big Bang 
event. The expected 4-volume of the cone cut at successive slices increases monotonically with depth. Fig. [To| :) 
shows that does not increase as fast as one might expect if each simplex on average had the same 4-volume (as in 
the unconstrained case): the number of nodes in each of the sub-volumes is #nodes(L) = 1, 5,15, 35, 70,126, 210, 330 
for L = 0,1,2, 3,4, 5, 6 , 7 respectively. Requiring causality thus changes the ’’shape” of the Universe even at strong 
coupling. Fig. 10 r) depicts the (unnormalized) distribution of the total 4-volume 14(7) of the lattice cone at /? = 1 


and 7 = —1.95. It resembles a F-distribution. 











































30 


Apart from an analytical continuation to real values of A > 0, the naive strong coupling limit of this model appears 
to be accessible numerically. Although devoid of curvature, this limit perhaps is not quite as trivial and uninteresting 
as might be presumed. 


X. SUMMARY AND DISCUSSION 


We considered the triangulation of causal Lorentz manifolds by a null lattice with hypercubic coordination. The 
nodes of this lattice represent space-time events. Spinors describe the light-like separation to neighboring events in 
the inertial system of the node at which they are based. Links between nodes are associated with SL{2, C) transport 
matrices. The basic 5'L(2, C)-invariant observables of this model are those of Eq.(|^). Despite the labeling by a 
lattice, this is a geometrical construction based on simplexes with four light-like edges and six spatial ones. It is similar 
to Regge ’s[8l|9] approach but retains fixed coordination and is based on an internal SL{2, C) symmetry. As shown in 
Sect|nj the triangulation of a causal manifold by light-like signals is physical and coordinate invariant. Macroscopically 
this method is used to triangulate events by radio signals cemitted from four GPS satellites. Conceptually the only 
difference here is that the emission events are themselves triangulated in the same fashion instead of being determined 
by atomic clocks on known orbits. 

We investigated two topologically different causal Lorentz manifolds: the universe with a single causally connected 
component and the duoverse with two. In some ways, the latter is aesthetically more appealing. Its lattice structure 
topologically is that of a body-centered hypercubic lattice whose two hypercubic sub-lattices are in dual position.The 
causally disconnected components of the duoverse are described by fields on these mutually dual sub-lattices which 
are joined by spatial links (on which SL(2, C) transport matrices live) only. Physical observables in the duoverse are 
constructed from closed lattice loops and the lattice Hilbert-Palatini action is given by a sum of elementary loops 
(see Eq.(43)). The duoverse could offer a novel accommodation for dark matter that does not interact causally 
with ordinary matter but does affect the observed local curvature. Depending on global boundary conditions, the 
two causally disconnected components of the duoverse in fact could be conjugate parts of a single non-orientable 
connected manifold and thus could lead to novel interpretations of the ’’arrow of time”, the absence of anti-matter 
in our part of the duoverse and other globally two-sided phenomena. However, a universe with a single causally 
connected component seems to be a viable, if less compelling, alternative. 

Ignoring global considerations, we first addressed local issues of the lattice theory. In Sect|V]necessary and sufficient 
conditions for a lattice configuration to represent a triangulated causal manifold were derived. Fo r nu ll lat-frames 
the local lattice ’’metric”, = |/^y(n)p, with (complex) amplitudes defined by Eq.(32al, is positive 

semi-definite, vanishing for fj, = n. The components of this symmetric matrix of spatial lengths satisfy the triangle 


inequalities of Eq.(35l with sides given by Eq.(36). These inequalities are equivalent to demanding that det(£^j,) < 0, 
which is necessary for a real local 4-volume. The consistency condition of Eq.(46) requires that the six spatial lengths 


on the backward light cone of a node satisfy the corresponding inequalities of Eq.(49), i.e. that the local volume of 


the corresponding backward null-simplex ch{n) is real as well. In the spinor formulation, this consistency requirement 
can be imposed by the TLT of Eq. ( [ 6 ^ . The corresponding TLT for null-latframes is obtained in Appendix 
Equivalently, these inequalities may be enforced by Heaviside distributions. 

The SL{2, C) structure group of this lattice theory is not compact and the invariant generating function diverges 
for any finite number of sites. In Sect VI the SL{2, C) structure group was partially localized to the compact SU{2) 


subgroup of spatial rotations. This is always possible and a local (equivariant) TFT was constructed that effects this 
partial gauge fixing. For simplicity and physical reasons, we mainly considered the constraints of Eq.( |78[ ). 

The integration measure of this lattice model is strongly constrained by its SL{2, C) structure group and local 
C7‘^(l)-invariance of the spinors. The measure for the SL{2, C) transport matrices is the (non-compact) Haar measure. 
The lattice integration measure for a spherical parametrization of the spinors was derived in Sect|VlTTj The invariant 
integration measure is determined up to a density p{V) that depends on local SL{2, C) invariants only, and in particular 
depends on the local 4-volume U(n). We used the gauge condition of Eq.([78| to localize this S'L(2, C)-invariant 
integration measure for the spinors. Although not necessary for numerical evaluation, we in addition completely 
localized the residual compact S'C/(2)-structure group. The spherical parametrization of the spinors in this case 
is given in Eq.( |108 ). This completely localized measure for the spinors in fact can be rewritten in terms of the 
manifestly SL{2, C) -invariant spatial lengths of Eq.(33) as in Eq.(lll). The restrictions on the length variables arise 
automatically in the change of variables. 

The apriori undetermined density p(V) of the integration measure allows for an invariant regularization of the 
short-distance behavior of the model that suppress es sm all local 4-volumes. Assuming the leading dependence of the 
density p(V) for V ^ 0 to b e power -like as in Eq. we argue that the lattice integration measure in the spinor 

parametrization of Eq.(108l is Eq.(llO). It depends on a parameter 7 . For a lattice with a large number N of sites, 
the dependence of 7 (A, N) on the cosmological coupling constant may be determined by demanding that the uni- (or 




















31 


duo-)verse have a given average total observable 4-volume. Critical values 7 * and A* are points of intersection of these 
curves as the number of sites iV —> 00 . 


Sect IX examined this model in the naive strong coupling limit, in which the cosmological term dominates and 
curvature corrections are ignored. The lattice action in this limit is ultra-local and depends on the local 4-volume 
only. Eq. ( [mt gives the local 4-volume in terms of the length variables. Since there is no dependence on transport 
matrices, the lattice model simplifies in this unphysical limit and the action and measure can be expressed in terms 
of the six invariant lengths at each site. We considered a conic lattice where the apex of the topologically 

hypercubic complex represents the Big Bang and the cone terminates at a spatial depth Tmax- We systematically 
found a number of additional non-compact surface symmetries that reflect the fact that there is some freedom in the 
placement of events on the boundary of the conic manifold: certain choices of events on the boundary can lead to 
different consistent triangulations of the same manifold and we therefore believe that these lattice symmetries are a 
consequence of residual diffeomorphisms on the boundary of the manifold. In Appendixj^this residual group of lattice 
symmetries is systematically derived. Since the group of residual symmetries is non-compact, it must be localized 
before attempting a numerical simulation. We find that this is achieved by imposing certain boundary conditions and 
adjusting the integration measure so that the expectation of invariant observables does not depend on this localization. 
In the naive strong coupling limit these observables include all partial lattice 4-volumes. The boundary conditions on 
spatial lengths of Eq.(E17) change the topology of the original hyperplane to depth L to that of a closed 3-manifold, 
much like the identification of the endpoints of a spatial segment to given depth L in the 1-1-1 dimensional results in 
a cone with topology R~^ x ^i. 

Implementing the additional boundary constraints, numerical simulations on small conic lattices with Lmax < 7 
and purely imaginary cosmological constant 4zA = 13 = 1 converge rather well in the naive strong coupling limit 
(see Fig. 10 They indicate that the critical exponent at which the 4-volume of the finite cone vanishes is 7 * ~ —2 
in the consistent model. This is to be expected since the consistency constraints are non-local, and 7 * = —2 when 
they are relaxed. The numerical results also show that the naive strong coupling limit of the consistent theory is not 
trivial: we computed the expected partial 4-volumes of the cone at depths 0 < L < Tmax = 7. Fig. 10:) shows that 


they increase much slower than the number of simplexes. When consistency conditions are imposed, the generating 
functional no longer factorizes into independent generating functions for each simplex and the simulations show that 
the correlations in fact are substantial. These initial numerical results are encouraging and perhaps warrant a more 
detailed investigation of the model. 

A number of issues have to be resolved before this lattice model becomes numerically accessible. Most importantly, 
the Lorentzian lattice weight of any configuration is a phase Lacking an analytic continuation, it is difficulfp^ 

to numerically evaluate expectations with complex generating functions. The oscillatory integrals in fact may be finite 
(as are Fresnel integrals) but are difficult to control and evaluate numerically. We so far only verified that the strong 
coupling limit of this model may be well defined for a purely imaginary cosmological constant. 

The complex weight of the Lorentzian model also makes it difficult to define an ’’average” proper time between 
two causally related nodes. Eq.(87) and Eq.( 88 ) are a reasonable definition of the proper time of a path and of 


the geodesic distance between two causally related active nodes of any particular causal configuration only. For an 
ensemble of such configurations, the proper time of a path between two particular nodes fluctuates and its average 
could even be a complex quantity. To address the general problem of defining distances in fluctuating geometries, 
one could consider conditioned amplitudes such as Eq.(90). In ordinary LGT’s this conditioning has no effect since 


distances between lattice nodes are known a priori on a space-time manifold that does not fluctuate. However, such 
a conceptual change in the definition of quantum amplitudes may be crucial in the construction of a quantum field 
theory based on ensembles of causal manifolds. 
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New approaches |454H8| . originally developed to accommodate a non-vanishing chemical potential in ordinary lattice gauge theories, may 
help overcome this problem. 
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Appendix A: Naturalness of the Cosmological Constant in Semi-Classical Quantum Gravity 


I argue semi-classically that the incredibly small value of the dimensionless cosmological constant A = ^ 

IQ-izo jg to the quantum fluctuations of the 4-volume V 4 of the observable universe. A similar argument is 

given in I3I1H1, bnt here I treat the cosmological constant as a coupling constant rather than the expectation of 
a fluctuating field. This coupling generally will itself depend on the epoch, i.e. on the 4-volume of the observable 
universe and on the energy density p of matter. The infrared free dimensionless coupling A = 327rA£p tends to vanish 
at large cosmic scales but A is assumed to remain large compared to ^pianckP cosmological term then dominates 
the evolution of the universe at later stages). Conversely, A could be large but dominated by the energy density in 
the early universe and lead to a strongly coupled theory that may be difficult to describe (semi-)classically. 

Let us therefore just note that fluctuations AV of the 4-volume of the observable universe proportional to the 
expected 4-volume {V), 

(AC)2 := (V^) - (V)^ = ijV , (Al) 


semi-classically imply a cosmological constant 


A = All = 32nAil 


Planck 


327 r(167r )£pianck 

ilVv 


(A2) 


Eq.(Al) follows if the total 4-volume of the universe is composed of a large number, N, of volumes Vi drawn 
independently from a single distribution whose variance is finite^ That the variance of the 4-volume is proportional 


to the 4-volume itself continues to hold under somewhat weaker assumptions. Much of the argument in is devoted 
to deriving Eq.( |Al[ ) from causal set theory [29ll4^ . Note that the present model comes close to a concrete realization 
of this idea: the total number of nearly independent 4-volumes in the strong coupling regime of Sect IX is A = 
at depth L where L is proportional to cosmic time at later stages of the evolution. 

Averaging over many oscillations of the phase ~ e-*AV4/(87r«p,,^„,,^) functional integral would lead to 

vanishing correlation functions}^ Fluctuations of the cosmological part of the action S thus cannot greatly exceed 
h = 27r/i. Neglecting curvature and matter contributions to the action, the cosmological term provides the main 
restriction on volume fluctuations. 

In epochs dominated by the cosmological part of the action, a semi-classical approach to quantum gravity therefore 
implies. 


A 


^’’’^Planck 


AE ~ 27r . 


(A3) 


With t he e stimate Eq.(AI) in Eq.(A3), one then arrives at Eq.(A2). 

In Eq.( A21 we introduced the Planck length £pianck = \/Gh/c^. It is the length scale associated with semi-classical 


gravitational phenomena and it would be unnatural if Ic <C ^pianck- For the purpose of a rough estimate it suffices to 
assume that ^ I'pianck, although it may be a few times larger. 

It remains to determine the average 4-volume of the observable universe. A somewhat naive estimate in terms of 
the radius r ~ 4.4 x I0 ^®to of the observable universe and its age cr ~ 1.3 x 10 ^®to yields 


A < (327r)16^y3£4j^„^i,/(4^r3cr) ~ 3 x IQ-^^o , (A4) 

which is about 1/5 of the value of A ~ 1.7 x 10 “^^® found experimentalljj^ and about twice (327r)/ty ~ 1.6 x 10 “^^° 
where tu = 8 x 10®° is the present expansion age of the universe in naturm units[3]. Considering the crudeness of the 
estimates, this is encouraging and an indication that gravity is quantized: A 10“^^° would imply highly suppressed 
quantum fluctuations and cast doubt on a semi-classical description of quantum gravity. A ^ 10“^^° on the other 
hand would imply an unnaturally large universe. 

Note that consistent with an infrared free coupling, quantum fluctuations of a 4-volume larger than the observable 
one ought to be larger and the cosmological constant correspondingly smaller. 

The cosmological constant often is considered an effective representation of the vacuum energy density. Its small 
value compared to typical (often UV-divergent) vacuum energies has been used as an argument for super-symmetric 


In this case // = {Av)^/v, where v is the mean and (An)^ is the variance of this distribution. 

The overall constant phase due to the average 4-volume cancels and is irrelevant in a semi-classical computation of amplitudes. 
The coupling A = AZp = 327rAtpj^^^j^ used here is about 327r ~ 100 times the cosmological constant in Planck units. 
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theories. However, the measured effective cosmological constant we are considering does not represent these vacuum 
energies. As with other coupling constants, only the renormalized cosmological constant can be measured. If quantum 
gravity is regularized, the measured effective cosmological constant is the difference of a bare cosmological constant 
and the ground state energy density of matter. It essentially determines the average 4-volume of the observable 
universe and has little to do with vacuum energies. 


Appendix B: First order Matter Actions 


The first order Hilbert-Palatini action of Eq.Q can be written in terms of forms. These diffeomorphism invariant 
objects retain their geometrical meaning when space-time is discretized. The coupling to matter also can be written 
in terms of forms and first order matter actions continue to depend polynomially on the co-frame 1-form (with the 
same sign ambiguity as the Hilbert-Palatini action). 

In curved space-time, the first order formulation of the invariant action for a gauge field with curvature two-form 
J-{A) = Fii;{A) dx^ A dx^ is. 


= ^ [ e“ A e*” A A e’^TvB^f B^f - 2TrH^‘^J^(A)]£. 

9 Jm 6 


abed 1 


(HI) 


where the trace is over internal degrees of freedom and Bab = —Bba are auxiliary 0-forms. A shift of the auxiliary 
fields leads to the equivalent and more familiar non-polynomial action. 


Sa = \ [ d^x det(e) Tr[F,fc(A)F*'=(A) - , with F*'=(A) = g^^g’^^FUA) , 

9 Jm 


(B2) 


if the co-frame is invertible. The first order action of Eq.(Bl) thus differs from the more conventional quadratic action 


by a local correction to the cosmological constant that depends on the auxiliary fields. In the lattice formulation the 
auxiliary fields could be integrated out. This would change the local density of the integration measure in Sect |VHr^ 
to by a factor proportional to the local 4-volume V~^, p(^) y~^p(y)- Correlators that do not depend on the 

auxiliary fields thus remain unaltered and the additional term does not change the dynamics and can be absorbed in 
a change of the local integration measu re. Note that Eq.( |B2[ ) again is prop ortional to det(e) rather than |det(e)|. 
Contrary to the quadratic action of Eq.(B2), the first order form in Eq.(Bl| is polynomial in all fields and depends 
on the co-frame 1-form only. 

There are some notable similarities between Eq.(Bl) and Eq.Q. Eq.( |Bl[ ) in fact reduces to the Hilbert-Palatini 
action for Tt:B°‘^F{A) = This would be the case for an so(3, l)-curvature form in the fundamental 

representation and an auxiliary field {B°‘^)cd = ~ ^d^c) 9^ = 2A/3. The action is stationary for this 


particular auxiliary field only if the curvature, A e° = 0 vanishes. The classical Hilbert-Palatini 

action thus can be interpreted as an so(3,1) gauge theory p^H22] - with pure gauge configurations of so(4,1) as classical 
solutions. Diffeomorphism invariance in this case is a property of the classical solutions only [21). 

The first order action for a scalar field in general coordinates is of the form, 


S^= [ e^ ATiilB^Bae^^ - B‘^d4>]eabcd 

J M 


(B3) 


The auxiliary fields again may be integrated out in favor of changing the density of the integration measure of 
Sect|VHIB|by p{V) V-^p{V). 


The invariant action of a spinor field may also be written in terms of the co-frame 1-form without explicit reference 
to frames. 


Sp = ( e“ A A e'^ A cr'^'D'tjj Sabed = [ d'^x det(e) ip ■ (B4) 

J M Jm 

Here '0 = '(/i a°''Dk'^P dx^ is a 1-form and {cr“, a = 1,..., 4} is a normalized basis of anti-hermitian 2x2 matrices. 

Note that the spinor action is polynomial only if the so(3, l)-connection of the covariant derivative 'Dip = dip+u}°'^aab'tp 
is assumed to not depend on the (co-)frame. 

Without cosmological term, the Jacobian of the invariant volume may be absorbed in the Hamiltonian formulation [5] 
of the Hilbert-Palatini action written in terms of frames. This is not possible when a cosmological term or matter 
interactions are included. In terms of frames, the action of the interacting quantum theory in this case is inherently 
non-polynomial. A first order formulation of the action using co-frames on the other hand remains polynomial when 
coupled to matter. 
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Another reason for quantizing co-frames rather than frames is that the length scale Ip can be absorbed in e^/Zp. 
The canonical mass dimensions (with h = c = 1) of all fields and couplings in this case is non-negative in d = 4 
dimensions. 


held or coupling 


. .ab 

A 

Afc 

^ab 

9 

0 


^l)1p 

can. mass dim. 

1 

1 

0 

1 

0 

0 

0 

0 

0 


Table 1: Canonical dimensions of helds and couplings in d = 4 dimensions 


The first order actions depend on the dimensionality d of space-time in two ways. The totally anti-symmetric 
Levi-Civita symbol e{abcd) with d arguments plays a prominent role in Eqs. ( l] |B4 ) and the number of co-frames in 
the monomials is determined by the order of the form they couple to [all actions may be written as inner products of a 
p-form constructed from the co-frames and a d — p form that does not depend on them]. The number of local d-forms 
one can construct with the given field content to some extent determines the form of the local action, although not 
as uniquely as desired since (j>, ipili and the auxiliary 0-forms are dimensionless. However, the local Hilbert-Palatini 
action of Eq.Q is unique in d = 4 dimensions. The other three local 4-forms one can construct with this held content, 
e“ A A i?ab(w), A and A R‘^‘^{uj)eabcd, are topological densities for the Nieh-Yan, Pontyagin 

and Euler characteristics of the manifold |35j that do not alter the continuum dynamics. 

In four dimensional space-time, vertices couple up to 6 helds in Eqs. (B1) (B3) and (B4) but at most 4 in Eq.([T|. 
In d > 3 dimensions, the hrst-order actions all lack quadratic terms. The lowest vertices in d = 4 couple 3 helds 
and there is no “natural” S'L(2, C)-invariant background about which a consistent perturbative expansion can be 
performed. 


Appendix C: Reconstruction of Backward Rays 

We here explicitly reconstruct four light-like vectors = (fp, —r^) with = jfpj > 0 on the (past) light cone of 
a node whose inner products. 


-2E„ 


Ejj — 2r^ ' 


71 


for p < i/ , 


(Cl) 


are known. The spatial lengths >0, l<p<i/<4} are assumed to satisfy the inequality Eq.(50). 


This reconstruction determines the local geometry of a lattice conhguration and is found to be unique up to (local) 
Lorentz transformations. Note that in the universe the lengths are the invariant spatial distances between two 
active diagonal nodes of a lattice plaquette whereas they give the displacement between two passive or virtual nodes 
in the duoverse. 

As explained in SectjVTjwe may choose a local inertial system in which. 




(C2) 


Squaring Eq.(C2) and using Eq.(Cl| to eliminate the scalar products, one finds that, 





(C3) 


where r denotes the average temporal separation of the four events from the 
minimal in the inertial system of Eq.(C2). Using that, 

4t(t^ +Ty-Tp- Ta) = 1 % - for {p, V, p, tj} = 
one obtains the (rotation invariant) temporal separations to the four events, 


node. Note that by construction r is 
{1,2,3,4}, (C4) 



(C5) 
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The condition of Eq. (C2 1 is invariant under rotation of the inertial system and we may choose one in which (in 
spherical coordinates), 


Ti = Ti(sin 0 i cos sin 9i sin <pi, cos 9i) 
r 2 = T 2 (sin 02 cos sin 6*2 sin (^ 2 , cos 6 * 2 ) 

= r 3 (sin 03 ,O, COS 03 ) 

fq = T4(0,0, 1) , 


(C 6 ) 


with 0 < 0j < TT and 0 < tpi < 27r. In this inertial system, the three angles 0,; a re uniquely determined in terms of the 
lengths £ 14,624 and £34 and the temporal components obtained in Eq.(C5l, 


sin( 0 j/ 2 ) = 


£i4 


2^^ 

The remaining angles ^pi and Lp 2 are found by noting that, 


< 1, i = l,2,3. 


(C7) 


cos ipi = 

cos ip 2 = 


13 


2 riT 3 (l — cos 01 cos 03 ) — £' 
2 rir 3 sin 01 sin 03 
2 t 2 T 3(1 — cos 02 cos 03 ) — 
2 r 2 T 3 sin 02 sin 03 


(C 8 ) 


Although the right hand sides are known, Eq.( |C 8 [ ) does not determine the angles 1^1 ^2 uniquely since cos(^ = 
cos(27r — if). One thus in general has four distinct solutions to Eq.(C 8 |. The choice is narrowed to (in general) two 
possibilities by requiring that, 


£^2 = 2 riT 2 (l — cos 01 cos 02 — sin 01 sin 02 cos((pi — P 2 )) ■ 


(C9) 


The reconstruction becomes unique for an orientable manifold with positive volume element, since the 4-volume 
spanned by the lat-frame is positive only if, 


cos 01 sin 02 sin 03 sin p 2 — cos 02 sin 03 sin 0i sin px + cos 03 sin0i sin 02 sin((/j 2 — <Pi) > 0 . 


(CIO) 


Orientability thus selects between the two remaining choices for pi^ 2 - 

We here assumed that the solutions to Eq.(C7l, Eq.(|C 8 |) and Eq.(C9) are consistent, i.e. that all angles are real. 


As argued in Sect ^ this reconstruction is possible if the six spatial lengths {£^ 1 / > 0, p < v} satisfy the inequality 
of Eq. (50). Note that we implicitly require this inequality in the reconstruction when the 4-volume is assumed to be 
real and positive. 


Appendix D: Consistency in terms of lat-frames 


We here construct a local TLT that imposes the constraints of Eq.(46) and leaves the additional iS'L(2, C) free. 
This equivariant TLT [SH] is based on the nilpotent BRST transformation. 


sEl{n) = c“(n) +a;“f,(n)£;|;(n) 

■sc“(n) = <()“f,(n).E|((n) +a;“6(n)c|((n) 

= w“c(n)<(>^b(n) - 

sa;“t(n) = - (j)%{n) 


sc^^n) = 
sg^^n) = 0 
sb^^''(n) 
sCo^^-'in) = 0 


and the BRST-exact action, 

Stlt = s X! [c'"''(n)(A;,(n - g. - v) ■ E^[n - g - v) - E^{n) ■ E^{n)) + 26 '"'^(n)E^(n) • c^{n)] 

= X! - g-y)- E^{-a. - g-v)- .E^(n) • E^{n)] + 

+ 2[c'"''(n) + w'"''(n)]E^(n) • c^{n) + 26'"‘'(n)[E^(n) • (j){n) ■ E„{n) + c^(n) • c^{n)] . 


(Dl) 


(D2) 
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Here c“(n) are the ghosts associated with the backward lat-frames E^(n). The ghosts a;““ = —w™ generate infinites¬ 
imal local SO~^{3, 1) transformations of the fields in Eq.(Dl). are the commuting topological ghosts [5^ 

that maintain the nilpotency of s. The remaining fields impose constraints and transform as doublets under this 
BRST. is a symmetric bosonic Lagrange multiplier field enforcing the constraints (46). c^^(n) = c‘'^(n) are 

the corresponding 10 anti-ghosts at each site. The are 6 bosonic fields with ghost number —2 that 

constrain the 16 ghosts c“(n). The a;^'^(n) = are their BRST-variations (with ghost number —1). The 

TLT-action of Eq.(D2| by construction is S'O’*'(3,1) gauge invariant and does not depend on the ghost that 
generates proper Lorentz transformations. 

Changing variables of the topological ghost to = s{E^ ■ c^) = E^(j)abE^ -I- • c^, the bosonic integrals over 

and 4’^'' lead^to a factor (det£’(n))“^ at each lattice site. Combining symmetric and antisymmetric anti-ghosts to 
= c^^n) -|-d;^^(n), there are altogether 16 anti-ghosts per site and a corresponding number of ghost fields 
c“. The integration over ghosts and anti-ghosts yields a factor (detE“(n))^. In any space-time dimension, this TLT 
imposes altogether d{d + l)/2 local constraints and leads to a factor detE“(n) in the measure at each site. Note that 
this determinant arises due to the quadratic dependence of the constraint on E)) - we could have deduced this factor 
in the measure by a Faddeev-Popov like procedure. For an orientable manifold one must choose the orientation of the 
Efj^ so that detE“(n) > 0 at every site. 

In the TLT one sim ply drops (or saturates) the integral over the ghosts a;)((n). Equivalently one may eliminate 
these fields in Eq.(Dl I from the outset and define an equivariant BRST transformation Se that is not nil-potent. 


SeL^;)(n) = c“(n) 
Sec“(n) = (t)\{x))El{n) 

s<?i>“&(n) = 0 


Seg^''{'o) = 0 

Seb^^''{n) = ^'^‘^(n) = -d;'^^(n) 

SedJ^''{n) = 0 . (D3) 


Sg is readily seen to generate a local SO'^{3, 1) transformation with parameters (n). 

Upon integration of the (j)^v,b^'',C^'' and c“-fields at each site we thus formally arrive at the local effective TLT 
action. 


9 '"'(n){Ef,{n - g-v)- E^{n - ^i-v)- E^(n) ■ E^(n)) + c^(n)U“(n)ca(n) 


(D4) 


where the integral over the new 4-component ghosts Ca(n) and anti-ghosts c^(n) yields the factor detU“(n) in the 
measure at each site. It is interesting to note that integrating some of the ghosts and Lagrange multiplier fields of 
the TLT action in Eq.( |D2[ ) led to the effective local action of Eq.(D4) whose eBRST-symmetry is no longer manifest. 
Integration of the (anti-)ghosts c^(n), Ca(n) and of the Lagrange multiplier g'^'^(n) gives the TLT partition function, 

Ztlt[E] oc /']^det.E(n)0[detU(n)] ]^d^U^(n) 0[-U))(n)] 5{E^{n- g,- v) ■ E^{n - g - v) - E^{n) ■ E^{n)) 

n fi 

oc ]jT(a(n),6(n),c(n)) ( f dP d{E^{n - g - v) ■ E^{n - g - v) - PU 

n / 




(D5) 


where a(n), 5(n), c(n) are given in terms of £‘^{n) = —2E^{n — g — v) ■ £'^(n — g — v) = £'^{n — g — ly) > 0 as in 


Eq.([56| and the distribution T(a,b,c), defined in Eq.(73l, imposes the non-local inequalities of Eq.(|49[). To arrive 


at the second line of Eq.(D5) we fixed the S'L(2, C)-invariance of the backward lat-frames as was previously done 
for the forward lat-frames. We then changed integration variables from backward lat-frames to the corresponding 
5'L(2, C)-invariant lengths. The resulting integration measure for the lengths is given by the analog of Eq.(lll) 
(with 7 = 0). The detU factor cancels against the inverse local volume V~^ of Eq.(lll). Due to the ^-distributions. 


the restrictions on the domain of the backward lat-frames become the non-local inequality constraints on the forward 
latframes imposed by T(a(n), 6(n), c(n)). 


The inequalities of Eq.(49l thus may be imposed by a TLT with the local action of Eq.(D4). 


29 


In d space-time dimensions, the Jacobian for this change of variables is det M oc (detE)"* ^ where M is the d(d — l)/2 X d(d — l)/2 
matrix — E“(n)E)j(n) constructed from the d X d-matrix E(n) at each site. 
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Appendix E: Boundary Conditions and Symmetries of a Conic Lattice 


1. Naive strong coupling without consistency constraints 


Ignoring the consistency constraints, the integration measure of Eq.(lll) at each site n may be written, 



(_5)(7-l)/20(_g) 


(El) 


where g is a a quaratic form in , 


g = det = a'^ + b'^ + — 2a^b^ — 2b^(? — 2(?a^ = (a^, 6 ^, (?) ^ ^ ) 


(E2) 


with a, 6 , c given in terms of the lengths as in Eq.( 112 1 . g depends only on the squares of spatial separations 

and the SL(2,R) isometry that leaves 5 < 0 invariant is generated by the real matrices. 


= 


^ ( 11 - 1 ) 

Vo 0 - 2 / 


rj-yZ _ 


1 


( 1 - 3-1 \ 
\2 -2 0 / 


2v^ V 2 -2 0 


_ 


1 


f-? 


0 1 -1 
0 1 


^3 V 1-1 0 


(E3) 


They commute to: [T^,T^]=T\ [T^,T^]=T^, [T\T^] = -T^. 

The action of this SL(2,R) on the 5-dimensional space, 

Q = {(fi 4 ,^ 241 '^ 34 j^ 23 i^i 3 i ^ 12)1 with det = — 1} , (E4) 


is induced by (for instance) considering Q as a fibre bundle over the 3-dimensional base B = {v := ii 3 >^i 2 )}> with 

2-dimensional fibers Fy := {/ = (^ 14 , 1 ^ 24 :^ 34 );5«(/) = 1} where 5 „(/) is the quadratic form of Eq.(E2) restricted to 
a base point, 


C _p4 p2 p2 p2 p2 \ 

^23 ^ 23'^13 ^ 23'^12 \ 

f-lAi (E5) 

p2 p2 p2 pi p4 I 

•‘- 12^23 •‘- 12^13 ^12 / 

The SL(2,R) acts on the fibers Fy of this bundle with generators, 

f«(i;) = diag(4-3^^^3^^^2') cliag(4,4,4) , (E6) 

and the corresponding variation of spatial separations is, 

= ff''{v)£l 4 with = 0 • (E7) 

To geometrically interpret the action of this SL(2,R), note that it changes only three lengths with a common vertex of 
the spatial tetrahedral face in a manner that preserves the 4-volume of the light-like 5-simplex. We chose a particular 
SL(2,R) subgroup of all 4-volume preserving transformations that leaves v := , ^\ 2 {'a)) invariant and 

acts on the three lengths / with common 4*^ vertex. We could just as well have considered isomorphic SL(2,R) groups 
that act on another vertex and do not change the lengths of the corresponding triangular base of the tetrahedron. 
We label these SL(2,R) groups by the vertex they change, as They are related because Q is invariant under 

3 continuous dilation symmetries that do not transform a,b or c. These dilations {D'',i = 1,2,3} commute with the 
and transform the lengths as. 


£14 




Hi4 , t 




23 


e “^^23 ; ^24 




H 24 , i 




13 


e-“=ri3 


£13 ; 44 e“434 , £12 ^ e-“4 


12 


(E8) 


It is readily verihed that any transformation is the composition of an transformation with two dilations 

D£ and where {i,j,k} S {1,2,3}. 

The volume preserving group of transformations of a null-simplex ch{n) therefore is isomorphic to 

g = X X X , (E9) 

and the remaining symmetry of the unconstrained lattice integration measure with N sites is . 
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The expectation of quantities invariant under can be evaluated by localizing the invariant measure of Eq.(El |. 
One for instance can demand at every node that, 


— hP' — and (P — 

a — U — C dlia — ^23, ^24 ~ ^135 "^34 “ ^12 5 


which is the same as requiring that the tetrahedron of the null-simplex be regular [50] with, 

•\/3£^,,(n) = y/—g(n) = E(n), for all 1 < /x < i/ < 4 . 


(ElO) 


(Ell) 


The localization conditions of Eq.(ElO) ca n be (uniquely) imposed. The last 3 conditions of Eq.(ElO) evidently are 
reached by the dilation symmetries of Eq.(E8). To show that = (? can always be satisfied by an SL(2,R) 

transformation, consider the Morse function, 


M{h) = a^{h) + b^{h) -b c^(h) > 0 


(E12) 


as a function of h gSL(2,R)/SO(2). M{h) is invariant under the compact SO(2) subgroup of SL(2,R). The critical 
points of M.{h) on SL(2,R)/SO(2) satisfy, 

0 = (a^ -b 6^ - 2c^){a^ + b^ + c^) = {a^ - b^){a^ + b^ + c^) . (E13) 


For -b 6^ -b > 0, M.{h) is minimal when a? = b^ = (? and the SL(2,R) transformation that reaches this 
minimum is unique up to (compact) SO{2) rotations. Invariants of Q thus are functions of the local 4-volume 


V = ~'\J~ = y/~9 only- 111® purpose of computing the expectation of an invariant, I{V), one thus may 

localize the 6-dimensional integral of Eq.(El) to a one-dimensional one over the 4-volume V of each null-simplex. 



(_g)(7-i)/2e(_^)j(y) 


dv V'^+^I{V) 


(E14) 


Without the consistency constraints imposed by T(a, &,c) in 
coupling limit would localize to. 


Eq.(114|, the generating function in the naive strong 


which is 


Eq.(116|. 


Zsc[A; 7] ot 

neA 


dV{n) E(n)^+^e 


7+1 —4iAy(n 


y-Af(7+2) 


(E15) 


2. Consistent naive strong coupling limit 


The factors T(a,6,c) = 0(—det (n)) in Eq.(114) enforce consistency constraints of the complex which ensure 


that every configuration represents a triangulated causal manifold. These constraints destroy many of the symmetries 
of the measure we have in their absence and the generating function no longer decomposes into independent factors for 
each node. Consistency thus imposes some correlations that prevent us from analytically calculating the generating 
function. In the end we have to resort to numerical simulations in the consistent naive strong coupling limit. 

However, a non-compact subgroup of symmetries S C survives the consistency constraints. The generating 
function Zgc diverges and an efficient numerical simulation of the consistent lattice model at strong coupling demands 
that these non-compact symmetries be localized. 

S is systematically found by considering the action of subgroups of Q on null-simplexes c/i(n € dC) on the boundary 
of the hypercubic cone C and extending their action in a manner that preserves the 4-volumes of all forward and 
backward null-simplexes of the lattice. We determine the full set of symmetries in this way because any spatial length 
£^iy(n) between two nodes, n -b /x and n -b ix, on a (cubic) face of the boundary, dC, is not constrained by consistency. 
The finite hypercubic cone C here is the collection of nodes. 


C = {n = ^ 




S N with L = 


4 

+=1 


/x=l,...,4} . 


(E16) 


At least one of the vanishes for nodes n G dC on the boundary of the cone C . There are four generic types of 
nodes to consider on dC: 
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1 ) 0 = ni = 712 = ns = n4. 

The apex of this null-simplex denotes the Big-Bang event. The diagonal links [n -1- /x, n -f z/] of c/i( 0 , 0 , 0 , 0 ) are 
all on faces of dC and their spatial lengths therefore are not constrained by consistency. The full group ^ is a 
symmetry of ch{0, 0,0,0). 

2 ) 0 = ni = 712 = ns, 774 > 0. 

This 1 -dimensional subspace of nodes is the 4*^ spine of the cone C. The diagonal links [n -|- /r, n -|- z/] of spinal 
null-simplexes are all on faces of dC and their spatial lengths are not constrained by consistency. A spinal 
simplex thus is invariant under the full group Q. 

3 ) 0 = Til = 772 ,713 > 0 ,714 > 0 . 

This is a 2 -dimensional subspace of nodes on the surface of C. All diagonal links except [n -|- Ai, n -|- A2] are 
on a face of the tetrahedron and therefore are not constrained by consistency. The corresponding length £12 (n) 
is invariant under x x C G. Note that and each act on a pair of lengths on faces of dC. 

The action of rescales the lengths £34(11) (which is on a face) and £i2(n) = £12(11 -t- Ai -|- A2) inversely. To 
preserve the volume of the (adjacent) backward null-simplex ch{n + Ai -|- A2), this dilation also has to rescale 
£34(n+Ai-|-A2) = £34(11-1-Aid-A2 —A3 —A4). It therefore also has to rescale £i2(n-f-A1-I-A2 —A3 —A4), etc. 
As discussed below, this continues until the length £12(11') between two nodes on another face is rescaled. The 
length £12(11') scales inversely to £34(11) on the original face. The action of on c/i(n) thus can be extended 
to a non-local symmetry that involves a string of cells. 

4) 77i = 0 ,772 > 0 ,773 > 0 ,774 > 0 . 

The diagonal links [n-|-A2, n-f A3], [n-f A2, n-f A4] and [n-f-A3, n-f A4] of length £23(11), £24(11) and £34(11) 
are on faces of dC but do not have a common vertex. The action of the on this null-simplex can be 

extended to a symmetry. Although it changes the lengths of £12(11),£13(11) and £14(11) which this cell shares 
with backward null-simplexes, the transformation may be followed by appropriate dilations that restore 

these lengths and instead rescale the separations £23(11),£24(n),£34(11) between nodes on faces. As in case 3 ), 
the action of D^,D^ and dilations on this cell can be extended to non-local symmetries that also rescale 
lengths of neighboring cells. 


Note that these four cases actually represent classes of surface nodes with four, three, two or one of the 77 ^ vanishing. 
They correspond to nodes at the intersection of 4, 3, 2 or 1 face(s) of the hypercubic cone. 

The remnant S C G^ symmetry group that survives the consistency constraints thus includes a set of non-local 
dilation sym metri es Dn\ associated with surface nodes n G dC of the cone. Their generators {Sn\i = 1,2,3} are 
given in Eg.(115 ). For {i,j, k} G { 1 , 2, 3}, Sn'^ generates a symmetry that scales £74(11') oppositely to £jfe(n'), just as 
U® in Eq.( |E 8 | ) does, but it scales these lengths on a one-dimensional string of lattice sites n' separated from n by an 
integer multiple 0 < p < q = min[ 77 i, 774 ] of Aj -|- A/j — A 7 — A 4 := These sites are all on the same depth L as 
n and in cases 1 ) and 2), q = 0 for all three I?^®\ whereas g > 0 for one dilation in case 3) and all dilations in case 
4). Note that there are exactly half as many independent dilations as there are diagonal links on the surface of the 
hypercubic cone. 


An SL(2,R) symmetry associated with each surface node also survives. T. Jacobson suggested [3^ that the residual 
group S of 4-volume preserving symmetries on the conic null-lattice is due to diffeomorphisms that relate triangulations 
of one and the same manifold. Even if interior nodes of the triangulation are determined by the intersection of light 
cones and form a null-lattice, some choice in the triangulation of nodes on the 3-dimensional light-like surface of a 
conic manifold remains. The residual group S relates equivalent null-lattice complexes with the same local 4-volumes. 
The last restriction implies that the symmetry group of equivalent null lattice configurations may be larger than 
S. However, for Im(A) < 0 the 4-volume of a null-simplex effectively is bounded and it suffices to consider only 
equivalence classes of configurations under the non-compact group S to define a finite generating function in the naive 
strong coupling limit. 

The residual scaling symmetries can be localized by imposing the boundary conditions. 


eUn)=i%{n') 


(E17) 


where the diagonal links [n -|- A 
n + min [777, 774 ]A(®\ with A^®^ given in Eq.(|ll5|. 


+ A 4 ] and [n' -|- Aj,n' -|- A*,] are on faces of the hypercubic cone and n' = 


The residual scaling symmetries allow the boundary condition of Eq.(E17) to always be satisfied. By viewing 
Eq.(E17) as a (unique) localization of the residual scaling symmetries, one can ensure that the expectation of invariants 
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formally does not depend on the boundary condition of Eq.(E17l by inserting, 


1= Yl ( 4 / - e 

[n+A,,n+A4]G9C ^ ^ 


(E18) 


in the generating function Zsc of Eq.(114|. Invariance of the measure (and of the strong coupling action and of 


observables), implies that the divergent volume of the residual scaling groups factorizes and may be absorbed in the 
normalization of the generating function. The product of integrals in Eq.(E18) extends over all sites n for which the 
diagonal link [n + , n + A 4 ] is on a face of, dC. 


Upon dropping the (divergent) volume of the residual scaling groups, one effectively has inserted, 

n (^?4(n)'5(^ffe(n + min[ni,n4]AW) - £-4(11)) 

[n+Ai ,n+A4]G5C 


(E19) 


in the integration measure of Eq.(114). 


We thus have imposed the boundary condition of Eq.(E17) while correcting the integration measure so that the 
expectation of invariant observables remain unaffected. Observables that are invariant under the scaling group in 
general are not local. They construction include any partial 4-volume of the cone and in particular the 4-volume of a 
slice of the lattice at a given depth L of the cone. 

As in Sect El the infinite volume of the SL{2,R) symmetry at each surface node also has t o be factorized and 
absorbed in the normalization before attempting a numerical simulation. As shown in Sect |E l[ SL(2,R) invariance 
can be used to achieve that, 


£i4(n)£23(n) = £24(n)£i3(n) = £34(n)£i2(n) 


(E20) 


for any null-simplex ch{n) with n S dC. 

As with the scaling symmetries, these conditions can be enforced without affecting the expectation of invariants by 
viewing Eq.(E20) as an equivariant localization of the coset SL{2, R)/SO{2). It can be shown that inserting. 


n (44<5 (^m^23 - - iU 

neac 


24^13 


34^12 


))„. 


(E21) 


in the integration measure of Eq.(114) achieves this. 


We localized the finite lattice model by inserting Eq.(E21| and Eq.(E19l in the integration measure of Eq.(114) 


and used the ^-distributions to evaluate the integrals of certain lengths of boundary nodes. The remaining lattice 
integral of the naive but now consistent strong coupling limit is finite and exhibits the correct scaling property: within 
numerical accuracy, the total expected lattice 4-volume is finite for 7 > 7 * = —2 for any cosmological constant with 
Im(A) < 0 and vanishes in the limit 7 —>■ 7 * = 2. Numerical results for several invariants on a small lattice are 
presented in Fig. [T^of Sect|IX[ 
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